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executive  summary 


I'he  Vortex  Advisory  System  (VAS)  was  developed  as  a means  for 
decreasing  costly  congestion  delays  in  aircraft  approach  queues  at 
the  major  airports.  The  system  evolved  from  the  analysis  of  the 
behavior  of  the  wake  vortices  from  over  50,000  landing  aircraft. 

It  was  determined  that  vortices  might  present  a threat  to  a follow- 
ing aircraft  only  during  specific  wind  conditions.  Most  of  the 
time,  however,  vortices  do  not  pose  a safety  problem.  The  VAS  in- 
dicates to  the  controllers  when  vortices  are  not  a problem;  during 
these  times,  the  interarrival  separation  standards  could  be 
reduced  to  3 nautical  miles  for  all  aircraft  regardless  of  leader 
or  follower  aircraft  type. 

The  vortex  behavior  data,  however,  were  collected  in  the  re- 
gion between  the  middle  marker  and  the  runway  threshold.  Over  82 
percent  of  the  landing  accidents  attributed  to  wake  vortices  oc- 
curred within  this  region.  The  VAS  advises  when  this  region  is 
clear  of  vortices.  Although  some  capacity  gains  might  be  realized 
if  the  reduced  separation  standards  were  used  only  in  this  region, 
the  utility  of  the  VAS  would  increase  if  the  protected  region  were 
extended  to  the  outer  marker. 

Applicable  vortex  behavior  data  did  not  exist;  therefore,  two 
programs  were  undertaken.  First,  the  subject  of  this  volume  of  the 
report,  an  analysis  was  done  to  calculate  the  relative  safety  of 
using  reduced  separations  from  the  outer  to  middle  markers  when 
the  VAS  indicated  that  such  reduced  separations  might  be  used  near 
the  runway.  Second,  the  subject  of  the  second  volume  of  the  report, 
data  were  collected  between  the  outer  and  middle  markers  to  verify 
the  adequacy  of  the  safety  analysis.  If  data  collection  alone 
were  to  be  used  to  establish  safety,  many  years  (and  many  data- 
collcction  sites)  would  be  required  to  collect  sufficient  data  to 
verify  the  safety  of  the  VAS.  However,  if  the  data  confirmed  the 
analytical  predictions  of  vortex  behavior,  the  analysis  herein 
would  demonstrate  the  safety  of  using  reduced  interarrival  spacings 
whenever  the  VAS  indicated  such  spacings  were  permitted.  Conversely, 
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ii  the  data  did  not  verify  the  analysis,  the  VAS  should  not  be  used 
to  reduce  spacings  between  the  outer  and  middle  markers. 

After  reviewing  vortex  - related  accidents  and  the  development 
and  implementation  of  the  VAS,  the  analytical  model  of  vortex  be- 
havior was  described.  State-of-the-art  knowledge  of  the  behavior 
of  vortices,  aircraft,  and  winds  aloft  was  used  to  formulate  the 
model.  At  all  times,  the  model  was  intentionally  caused  to  err  on 
the  conservative  side.  Thus,  the  risks  calculated  using  this  model 
were  proxy  estimates  of  the  real  stiuation 

The  risks  accepted  in  the  present  system  using  the  3-,  4-,  and 
5-nautical-mile  spacings  in  effect  near  the  outer  marker  were  cal- 
culated. Twelve  aircraft  were  considered  as  vortex-generating  or 
leading  aircraft,  and  the  same  12  were  considered  as  the  vortex- 
encountering  or  following  aircraft:  B-747,  L-1011,  PC-10,  B-707H, 

DC  - 811 , B-  707  , DC-8,  B-727  , DC-9,  B-737  , Learjet,  and  PA-28.  Accord- 
ing to  the  model,  at  the  current  separations,  the  risk  of  a hazard- 
ous-vortex encounter  is  zero  for  Large  and  Heavy  aircraft  following 
other  Large  and  Heavy  aircraft  and  nonzero  only  for  some  Small- 
category  aircraft  following  certain  Large  and  Heavy  aircraft.  The 
situations  represented  by  the  nonzero  risks  occur  frequently  in 
the  current  system  in  the  vicinity  of  the  outer  marker.  Since  these 
situations  occur  frequently  (particularly  in  VFR  conditions)  and 
with  nonzero  risk  of  a vortex  hazard,  these  risks  were  considered 
to  be  acceptable- -no  hazardous -vortex  encounters  or  accidents  were 
known  to  occur  for  these  situations. 

The  largest  of  the  nonzero  risks  was  defined  as  a baseline 
value.  If  the  chance  of  a hazardous -vortex  encounter  for  any 
1 eader/ fol lower  pair  were  found  to  be  less  than,  or  equal  to,  the 
baseline  value,  the  situation  would  be  considered  safe. 

Risks  to  be  expected  when  3-nautical-mile  spacings  would  be 
permitted  were  calculated.  Operational  guidelines  were  then 
formulated  to  maintain  risks  with  VAS-reduced  spacings  at,  or  below, 
the  baseline  value.  These  guidelines  are:  (1)  reduced  separations 
to  be  used  only  when  the  VAS  indicated  that  conditions  permitted 
such  separations,  and  (2)  precision  approaches  are  required  (i.e., 
no  short  finals  or  VOR/ loca 1 i zer  approaches).  Kith  these  guide- 
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lines,  and  the  verification  of  the  model  (Volume  II)  the  risks 
with  VAS- reduced  spacings  would  be  no  larger  than  in  the  current 
system  with  the  present  separation  standards. 

A brief  sensitivity  analysis  was  done  to  determine  if  the 
conclusions  of  the  safety  analysis  were  adversely  affected  by  pos- 
sible errors  in  the  models  and/or  parameters.  Although  the  details 
were  altered,  the  conclusions  concerning  the  safety  of  VAS-reduced 
spacings  were  not  changed. 

Other  evidence  was  investigated  to  look  for  support,  or  con- 
tradiction, of  the  analytic  results.  Wake  vortex  incident  reports 
analyzed  by  the  British  Civil  Aviation  Authority,  analysis  of  vor- 
tex-caused accidents,  and  the  current  spacings  actually  flown 
safely  in  VFR  all  lent  credence  to  the  results  of  the  calculations: 
VAS-reduced  spacings  should  be  safe. 

There  were  two  unresolved  issues  in  the  analvs is  - - the  model 
describing  cross  winds  aloft  is  unproven  and  its  derivation  re- 
quired somewhat  inconsistent  assumptions,  and  a possible  problem 
might  have  been  introduced  by  comparing  two  conservative  estimates 
of  the  probability  of  a hazardous-vortex  encounter.  Resolution  of 
the  former  issue  will  be  addressed  in  Volume  II  and  the  conse- 
quences (if  any)  of  the  latter  issue  in  Volume  III.  Assuming  favor- 
able resolution  of  these  two  issues,  the  analysis  indicates  the 
safety  of  VAS-reduced  spacings  between  the  outer  and  middle  markers 
lor  all  aircraft.  The  safety  of  VAS-reduced  spacings  between  the 
middle  marker  and  touchdown  was  established  by  the  study  of  the 
vortices  from  over  50,000  landing  aircraft. 

An  alternative  analysis  was  done  without  recourse  to  either  of 
the  two  unresolved  issues.  Assuming  that  50  percent  of  available 
roll-control  authority  was  acceptable  for  countering  the  rolling 
motion  imposed  by  a vortex  (higher  percentages  are  accepted,  or 
appear  to  be,  with  today's  separation  standards),  vortex  decay  would 
be  the  primary  mechanism  precluding  hazardous-vortex  encounters 
between  the  outer  and  middle  markers.  Vortex  decay  alone,  however, 


Vortex  decay  alone  was  shown  to  permit  VAS-reduced  spacings  for 
Heavy  and  Large  aircraft  following  Heavies. 
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PREFACE 


The  concept  of  a Vortex  Advisory  System  (VAS)  evolved  from  the 
analysis  of  an  accumulating  wealth  of  vortex-tracking  data.  Tens 
of  thousands  of  tracks  demonstrated  that  the  concept  was  viable; 
however,  all  the  data  were  collected  between  the  middle  marker  and 
the  runway  threshold.  Before  certifying  the  VAS  for  an  operational 
test,  the  Flight  Standards  Service  of  the  Federal  Aviation  Adminis- 
tration (FAA)  requested  that  the  region  between  the  middle  and  outer 
marker  be  examined  to  determine  its  relative  safety  vortexwise  for 
the  use  of  decreased  separations. 

Many  people  contributed  to  this  report.  The  mathematical  model 
was  developed  from  an  unpublished  working  paper  by  Alan  Bilanin  and 
Coleman  Donaldson  of  Aeronautical  Research  Associates  of  Princeton. 

It  is  a pleasure  to  acknowledge  their  work  as  well  as  the  able 
assistance  of  Berl  Winston  of  the  Transportation  Systems  Center 
(TSC)  who  programmed  the  model  and  helped  to  run  the  numerous  cases. 

I wish  especially  to  thank  R.  M.  Harris  of  Mitre  Corporation;  his 
critical  reading  of  an  early  draft  of  this  report  prompted  a com- 
plete revision  of  the  document.  He  derived  the  cross-wind  probability 
density  model  presented  in  Appendix  C and  made  many  detailed  sug- 
gestions for  clarifying  the  model.  Again  Sinha  of  Mitre  Corporation 
pointed  out  the  possible  problem  of  comparing  two  conservative  esti- 
mates (Section  5.2),  and  I wish  to  thank  him  for  his  detailed  com- 
ments which  are  reflected  in  Sections  5.2  and  6.1.  I thank  R. 

Craig  Goff  for  supplying  the  AIDS  data  which  are  discussed  in  Sec- 
tion C.4.  Finally,  I would  like  to  thank  Myron  Clark,  Joe 
Tymczyszyn,  Jr.,  Andrew  Haines,  and  Dave  Burnham  for  their  sug- 
gestions and  comments  on  the  final  draft  of  the  report. 
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1.  INTRODUCTION 


All  aircraft  gene  rate  trailing  wake  vortices  as  a consequence 
of  generating  lift;  however,  the  potential  danger  of  encountering 
these  wake  vortices  has  only  recently  become  apparent.  Aircraft 
wake  vortices  now  constitute  one  of  the  major  problems  confronting 
the  air  traffic  control  system. 

Before  1970,  landing  aircraft  were  required  to  maintain  at 
least  3-nautical-mile  separations  under  Instrument  Flight  Rule 
(IFR)  conditions.  The  separation  standard  was  based  primarily  on 
radar- operat ing  limits,  and  to  a lesser  extent,  on  runway- occupancy 
limitations.  There  were  no  separation  standards  imposed  because  of 
vortex  considerations. 

With  the  introduction  of  the  wide-body  jets  and  the  increasing 
number  of  aircraft  operations  at  the  major  airports,  the  wake  vor- 
tex problem  has  taken  on  increasing  importance.  The  vortices  from 
large  aircraft  can  present  a hazard  to  smaller  aircraft  which  inad- 
vertently encounter  the  vortices;  the  following  aircraft  can  be 
subjected  to  rolling  moments  which  exceed  the  roll  control  author- 
ity of  the  aircraft,  to  a dangerous  loss  of  altitude,  and  to  pos- 
sible structural  failure.  The  probability  of  a vortex  encounter  is 
greatest  in  the  terminal  area  where  Small,  Large,  and  Heavy  aircraft 
operate  on  the  same  flight  paths  in  close  proximity,  and  where 
recovery  from  an  upset  may  not  be  possible  because  of  the  low' 
aircraft  altitude. 

Accordingly,  the  solution  implemented  by  the  Federal  Aviation 
Administration  (FAA)  in  March  1970  was  to  increase  the  separation 
standards  behind  the  Heavy  jets  (a  Heavy  jet  has  a maximum  certifi- 
cated takeoff  weight  of  300,000  pounds  or  more)  evolving  by  1973  to 
4 nautical  miles  for  a following  Heavy  aircraft,  and  to  5 nau- 
tical miles  for  a following  non-Heavy  aircraft.  The  standards  were 
revised  in  November  1975  by  requiring  the  addition  of  an  extra  nau- 
tical-mile separation  at  threshold  for  following  aircraft  with  a 
maximum  certificated  takeoff  weight  of  less  than  12,500  pounds 
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behind  a Large  or  Heavy  aircraft.  These  increased  separations 
led  to  additional  delays,  and  decreased  the  capacity  and  efficiency 
of  the  airport  system. 


One  technique  that  has  been  developed  to  regain  some  of  the 
lost  capacity  is  the  use  of  a Vortex  Advisory  System  (VAS) . The 
system  indicates  to  controllers  when  the  separation  standards  can 
be  reduced  to  3 nautical  miles  regardless  of  the  leader  or  follower 
aircraft  type.  The  VAS  has  evolved  from  the  examination  of  tens  I 

of  thousands  of  vortex  tracks  and  the  correlation  of  vortex  behavior 
with  the  ambient  winds. 

All  the  vortex-tracking  data  have  been  recorded  in  the  middle 
marker  to  runway  threshold  region.  Over  82  percnet  of  the  landing 
accidents  attributed  to  wake  vortices  have  occurred  within  this 
region  (Ref.  1).  The  use  of  the  VAS  advises  when  this  region  is 
clear  of  vortices.  Although  some  capacity  gains  may  be  realized  if 
only  this  region  is  permitted  to  use  the  reduced  separation  stand- 
ards, the  utility  of  the  VAS  increases  if  the  protected  region 
can  be  extended  to  the  outer  marker. 

Since  applicable  data  did  not  exist,  two  approaches  were  under- 
taken. First,  the  subject  of  this  volume  of  the  report,  a proba- 
bility analysis  was  done  to  calculate  the  relative  safety  of  re- 
duced separations  out  to  the  outer  market  when  the  VAS  indicated 
that  reduced  separation  would  be  permitted  near  the  runway.  Such 
an  analysis  was  undertaken  as  it  was  expected  that  the  conditions 
that  permitted  reduced  separations  inside  the  middle  marker  would 
also  permit  reduced  separations  to  be  used  inside  the  outer  marker. 

The  second  approach,  the  subject  of  the  second  volume  of  this 
report,  concerns  the  gathering  of  appropriate  data  to  verify  the 
claims  of  the  probability  analysis. 

The  method  employed  in  the  probability  analysis  is  to  use  the 
best  data  available  which  describe  vortex  and  aircraft  behavior  in 
the  approach  region  to  calculate  the  chance  of  a vortex  encounter. 

At  all  times,  the  model  is  deliberately  caused  to  err  on  the  con- 
servative side.  To  give  meaning  to  small  probability  values,  the 
numbers  are  interpreted  by  comparing  with  a known  safe  situation. 
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In  general,  the  probabilities  calculated  in  this  report  have  been 
overestimated  in  absolute  terms,  but  the  relative  values  of  the 
calculated  probabilities  should  not  be  affected. 

Section  2 reviews  the  history  of  vortex  - related  accidents  to 
put  the  FAA  wake  vortex  program  into  proper  perspective.  The  most 
frequent  type  of  vortex- related  accident  involves  one  aircraft  on 
approach  following  another  aircraft  landing  on  the  same  runway. 

Most  of  these  accidents  have  occurred  within  the  middle  marker  and 
with  interarrival  separations  less  than  3 nautical  miles. 

Section  3 describes  the  development  and  embodiment  of  the  VAS 
with  particular  reference  to  the  Chicago  O'Hare  system.  The  key 
to  the  VAS  is  a wind  criterion  which  has  evolved  from  the  analysis 
of  over  50,000  aircraft  landings.  Whenever  the  surface  wind 
exceeds  a defined  wind  criterion,  IFR  interarrival  spacings  may  be 
safely  reduced  to  a uniform  3 nautical  miles;  whenever  the  surface 
wind  does  not  exceed  the  wind  criterion,  the  separations  remain 
unchanged  at  the  vortex-based  3,  4,  5,  or  6 nautical  miles.  The 
analysis  of  the  over  50,000  vortex  tracks  is  the  foundation  for 
the  safety  of  the  VAS  in  the  runway  threshold  to  middle-marker 
region . 

The  vortex  encounter  model  is  derived  in  Section  4.  State-of- 
the-art  knowledge  of  the  behavior  of  vortices,  aircraft,  and  winds 
is  used  to  construct  an  analytical  model  for  calculating  the  prob- 
ability of  a hazardous-vortex  encounter  between  the  outer  and 
middle  markers.  By  intent,  the  model  is  both  as  general  as  possible 
and  conservative.  The  development  of  the  model  mixes  rigorous 
mathematics  with  many  engineering  approximations.  Sufficient  data 
exist  to  permit  constructing  defensible  descriptions  of  vortex  and 
aircraft  behavior.  Imperfect  knowledge  of  the  correlation  between 
cross  winds  aloft  and  VAS-measured  winds  has  led  to  a very  mathe- 
matical, and  sometimes  inconsistent,  description  of  the  winds; 
limited  data  and  observations  seem  to  support  the  description,  but 
the  measurements  to  be  discussed  in  Volume  II  are  required  for 
verification. 

Section  5 contains  the  safety  analysis  for  the  outer  to  middle 
marker  region,  both  the  logic  and  the  results,  encompassing  all 
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aircraft  categories (Heavy , Large,  and  Small).  It  is  shown  that 
interarrival  separations  may  be  set  at  a uniform  3 nautical 
miles  in  I1:R  whenever  the  VAS  indicates  that  reduced  separations 
may  be  employed.  Aircraft,  however,  must  be  restricted  to  preci- 
sion straight- in  approaches  from  the  outer  marker.  A possible 
problem  is  indicated  concerning  the  comparison  of  conservative 
estimates;  the  analysis  of  this  problem  is  extensive  and  is  de- 
ferred to  Volume  III  along  with  a more  detailed  sensitivity  anal- 
ysis. Assuming  successful  resolution  of  the  problem  and  valida- 
tion of  the  wind  model,  the  analysis  in  Section  5 shows  why  3-nau- 
tical-mile  separations  can  be  safe  vortexwise  regardless  of  the 
category  of  the  leader  and  follower  aircraft. 

There  are  2 main  unresolved  issues- -the  cross-wind  aloft 
model  and  the  possible  problem  related  to  the  comparison  of  con- 
servative estimates.  In  Section  6,  the  scope  of  the  VAS  is  reduced 
to  apply  only  to  Large  and  Heavy  aircraft  landing  behind  Heavy  air- 
craft. Here,  the  safety  of  3-nautical-mile  separations  is  dis- 
cussed without  using  either  of  the  2 unresolved  issues. 

Sections  2,  3,  and  4 are  not  dependent  on  material  in  other 
sections.  Sections  5 and  6 refer  often  to  Section  4.  Depending 
on  the  interest  of  the  reader,  Sections  2 and  3 may  be  skipped. 

If  the  details  of  the  vortex  encounter  model  are  accepted;  Sections 
4.1.1,  4.1.2,  and  all  of  5 and  6 contain  the  material  necessary 
for  understanding  why  it  is  claimed  that  the  VAS  will  permit  safe 
3 -naut ical -mil e separations  between  the  outer  and  middle  markers. 

Section  7 summarizes  briefly  the  alternative  safety  anal- 
yses discussed  in  Sections  S and  6.  Both  analyses  use  the  same 
model  for  a vortex  hazard.  If  a maximum  roll-control  input  of  50 
percent  of  the  full  roll-control  authority  of  an  aircraft  to 
counter  a vortex  swirling  flow  is  acceptable  (assuming  the  worst 
case  of  an  axial  vortex  encounter  with  a very  slowly  decaying  vor- 
tex), then  3-nautical-mile  separations  between  outer  and  middle 
markers  may  be  used  for  Large  and  Heavy  aircraft  following  Heavies. 
To  extend  the  VAS  to  include  Small  aircraft,  two  unresolved 
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issues  have  been  introduced  into  the  analysis  (wind  model  and 
possible  problem  of  comparing  conservative  probability  estimates). 
If  favorably  resolved  (Volume  II  and  III),  then  the  VAS  can  be 
used  lor  all  aircraft  as  is  shown  in  Section  5. 


2,  WAKE  VORTEX  ACCIDENTS 


Of  the  approximately  45,000  aviation  accidents  that  occurred 
during  the  10-year  period  of  1964-1973  in  the  conterminous  United 
States,  wake  vortices  were  cited  by  the  National  Transportation 
Safety  Board  (NTSB)  as  a cause  or  factor  in  147  accidents.  Some  of 
the  147  accidents  may  not  be  genuinely  vortex-related;  analyses  sub- 
sequent to  the  issuance  of  the  NTSB  reports  question  the  role  of 
vortices  in  at  least  63  of  these  accidents.  Of  the  remaining  84 
probable  vortex- related  accidents,  27  resulted  in  fatalities. 
Reference  1 describes  a study  of  the  vortex- related  accidents,  and 
categorizes  these  accidents  by  the  flight  phases  of  the  vortex- 
generating aircraft  and  the  accident  aircraft,  their  relative  run- 
ways of  operation,  and  other  pertinent  factors. 

Examination  of  the  vortex- related  accidents  provides  historical 
perspective  to  the  wake  vortex  program.  The  discussion  below  fol- 
lows closely  the  summary  of  Ref.  1 published  in  Ref.  2.  The  statis- 
tics herein  update  Refs.  1 and  2 to  reflect  recent  findings  con- 
cerning two  of  the  accidents. 

2.1  OVERALL  STATISTICS 

During  the  period  of  1964-1973,  the  NTSB  classified  an  average 
of  15  accidents  per  year  as  being  vortex-related;  3 accidents  per 
year  resulted  in  fatalities  (Ref.  1).  Single- aircraft  accidents, 
by  comparison,  occurred  at  an  annual  rate  of  4510  with  540  fatal. 
Vortex-related  accidents  constituted  only  0.33  percent  of  all  single- 
aircraft accidents.  Twenty-seven  of  the  147  accidents  were  fatal 
and  resulted  in  55  fatalities  (Fig.  1).  There  were  88  landing 
accidents  (21  fatal  and  43  fatalities),  28  takeoff  accidents  (4 
fatal  and  9 fatalities),  and  31  in-flight  accidents  (2  fatal  and  3 
fatalities).  Twenty-six  of  the  31  in-flight  accidents  were  crop 
dusters  (1  fatal  and  1 fatality)  involved  in  agricultural  activities 
in  close  proximity  to  the  ground. 

Eliminating  the  apparently  vortex-unrelated  accidents  and  the 
in-flight  crop-duster  accidents  leaves  64  landing  accidents 
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(16  fatal  and  33  fatalities),  18  takeoff  accidents  (3  fatal  and  8 
fatalities),  and  2 in-flight  accidents  (1  fatal  and  2 fatalities). 

The  probable  vortex-related  accidents  are  84,  20  fatal  with  43 
fatalities  (Fig.  2). 

Approximately  five-sixths  of  the  landing  accidents  and  three- 
quarters  of  the  takeoff  accidents  occurred  at  controlled  airports. 

For  the  ten-year  period  of  1964-1973,  the  vortex-related  accident 
rate  at  towered  airports  was  less  than  1 per  3-million  landings  and 
less  than  1 per  10-million  takeoffs.  The  air  carrier  rate  was  less 
than  1 per  40-million  operations  for  both  landings  and  takeoffs. 

General  aviation  aircraft  weighing  less  than  12,500  pounds 
have  been  the  primary  victims  of  the  vortex  problem.  The  general 
aviation  aircraft  seems  to  be  almost  as  vulnerable  to  the  vortices 
of  aircraft  weighing  more  than  100,000  pounds  as  to  the  vortices 
of  aircraft  weighing  more  than  300,000  pounds  (Ref.  1). 

There  have  been  only  two-probable  vortex-related  air-carrier 
accidents:  a DC-9  landing  accident  (a  training  flight  behind  a DC- 
10)  at  Fort  Worth,  Texas;  and  a DHC-6  (Twin  Otter)  takeoff  accident 
at  Kennedy  International  Airport,  Jamaica,  New  York.  The  latter, 
however,  was  an  intersect  ion- takeoff  accident  involving  an  aircraft 
weighing  less  than  12,500  pounds. 

The  study  of  the  NTSB  data  base  indicated  that  all  the  probable 
vortex-related  landing  accidents  occurred  when  the  aircraft  separa- 
tions were  less  than  the  separation  standards  (Ref.  1) . The  only 
probable  vortex-related  landing  accident  that  occurred  under 
Instrument  Flight  Rules  (IFR)  was  an  in-flight  vortex  encounter  where 
both  aircraft  were  conducting  missed  approaches. 

2.2  LANDING/LANDING:  SAME  RUNWAY  ACCIDENTS 

The  most  frequent  type  of  vortex-related  accident  involved  an 
aircraft  on  approach  following  another  landing  aircraft  landing  on 
the  same  runway.  Fifty-seven  such  accidents  occurred  behind  con- 
ventional aircraft  (45  probable),  and  7 occurred  behind  landing 
helicopters  (all  probable).  This  category  thus  accounts  for  more 
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FIGURE  2.  BREAKDOWN  OF  PROBABLE  VORTEX-RELATED  ACCIDENTS  FOR 
THE  TEN-YEAR  PERIOD  1964-1973 


than  82  percent  of  the  vortex- related  landing  accidents,  and  about 
76  percent  of  all-probable  vortex  - related  accidents.  Since  the 
Vortex  Advisory  System  (VAS)  seeks  to  minimize  delays  by  decreasing 
landing/landing  separations  during  certain  wind  conditions,  the 
vortex-related  landing/landing  accidents  are  examined  in  greater 
detail.  Only  the  fixed-wing  aircraft-caused  accidents  are  consid- 
ered below. 

2.2.1  Encounter- Point  Statistics 

Table  1 presents  the  encounter-point  statistics.  For  the 
probable  vortex-related  landing  accidents,  33  percent  of  the  vortex 
encounters  took  place  after  the  victim  aircraft  had  crossed  the 
runway  threshold,  54  percent  between  the  runway  threshold  and  the 
middle -marker  region,  and  the  remaining  13  percent  occurred  more 
than  1 nautical  mile  from  the  runway  threshold.  Thus,  87  percent 
of  the  landing  accidents  occurred  within  the  region  of  the  middle 
marker . 

The  accidents  outside  the  region  of  the  middle  marker  seem 
to  be  rather  random  encounters  occurring  because  a small  aircraft 
is  approaching  low  and  with  poor  lateral  navigation.  In  this 
manner,  the  small  aircraft  manages  to  encounter  vortices  even  when 
the  vortices  have  been  blown  away  from  the  extended  runway  center- 
1 ine . 

Comparatively  few  accidents  resulted  from  vortex  encounters  at 
high  altitudes.  This  may  be  partially  explained  as  higher 
encounter  altitudes  provide  space  in  which  the  aircraft  can  recover. 

2.2.2  Aircraft  Types  Involved 

Figure  3 shows  a landing  weight  histogram  for  accident  air- 
craft involved  in  landing  accidents  behind  fixed-wing  aircraft 
landing  on  the  same  runway.  It  is  clear  that  small  general  aviation 
aircraft  is  the  primary  victim  of  the  vortex  problem.  About  85 
percent  of  the  accident  aircraft  weighed  less  than  4000  pounds, 
while  almost  97  percent  weighed  less  than  11,000  pounds. 
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FIGURE  3.  LANDING  WEIGHT  HISTOGRAM  FOR  ACCIDENT  AIRCRAFT  DURING  LANDING 
ACCIDENTS  BEHIND  LANDINGS  ON  THE  SAME  RUNWAY 


As  for  the  vortex-generating  aircraft  which  caused  these  acci- 
dents, the  largest  number  of  accidents  (12J  were  attributed  to 
B-707's,  followed  by  6 accidents  attributed  to  B-727's.  Figure  4 
indicates  the  weights  of  the  vortex-generating  aircraft.  More 
than  90  percent  of  the  accidents  have  been  caused  by  non-Heavy 
air-carrier  aircraft  weighing  less  than  300,000  pounds;  this  may 
not  be  significant  given  the  relatively  few  Heavy  aircraft  com- 
pared to  the  many  non-Heavy  air-carrier  aircraft  during  the  period 
of  this  analysis.  Some  accidents  have  even  been  attributed  to 
vortices  from  aircraft  weighing  less  than  25,000  pounds. 

2.2.3  Reported  Wind  Conditions 

Approximately  two- thirds  of  the  accidents  occurred  under  wind 
conditions  of  5 knots  or  less.  Some  accidents  in  the  NTSB  data 
base  did  occur  with  winds  higher  than  10  knots;  most  of  these  ac- 
cidents, however,  were  found  to  be  unrelated  to  the  vortex  problem. 

2.3  HISTORICAL  PERSPECTIVE 

The  preponderance  of  vortex-related  accidents  for  landing 
aircraft  occurred  between  the  middle-marker  location  and  the  run- 
way threshold.  Accordingly,  the  United  States  program  on  aircraft 
wake  vortices  (Ref.  2)  concentrated  on  the  middle  marker  to  thresh- 
old region  to  develop  a Wake  Vortex  Avoidance  System  (WVAS)  as  an 
integral  part  of  the  Upgraded  Third  Generation  Air  Traffic  Control 
System  (Ref.  3).  The  Vortex  Advisory  System  (VAS)  concept  evolved 
from  the  analysis  of  approximately  50,000  vortex  tracks,  and  is 
proposed  as  an  interim  measure  to  decrease  delays  at  the  major  hub 
terminals . 
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CD  NOT  BELIEVED  TO  BE  VORTEX  - RELATED 


II  PROBABLE  VORTEX -RELATED  ACCIDENTS 


NOTE:  "UNIDENTIFIED  .JETS"  NOT  INCLUDED  IN  THIS  HISTOGRAM 


FIGURE  4.  LANDING  WEIGHT  HISTOGRAM  FOR  VORTEX-GENERATING 
AIRCRAFT  DURING  LANDING/LANDING  ACCIDENTS  ON  THE  SAME  RUNWAY 
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3,  VORTEX  ADVISORY  SYSTEM 


Trailing  vortex  wakes  of  large  aircraft  can  pose  a hazard  to 
smaller  following  aircraft.  As  a result,  large  separation  dis- 
tances have  been  imposed  to  maintain  safety.  (Separation  distances, 
in  this  report,  refer  to  the  interarrival  spacings  between  air- 
craft imposed  during  final  approach  under  Instrument  Flight  Rules.) 
Under  the  sponsorship  of  the  FAA,  the  TSC  has  collected  an  abundant 
amount  of  information  about  vortices.  Analysis  of  the  extensive 
data  on  vortex  behavior  as  a function  of  meteorological  conditions 
has  indicated  that  there  are  wind  conditions  during  which  vortices 
do  not  pose  a threat  to  a following  aircraft  (Refs.  2,  4,  and  5J . 

3.1  VORTEX  DATA  COLLECTION 

After  limited  tests  had  been  conducted  at  Boston's  Logan 
International  Airport,  the  .John  F.  Kennedy  International  Airport 
became  the  first  extensive  test  site  for  the  study  of  vortex  be- 
havior. The  purpose  was  to  evaluate  two  vortex  - sens ing  systems, 
the  Ground  Wind  Vortex- Sens ing  System  (GIWSS)  and  the  Pulsed 
Acoustic  Vortex-Sensing  System  (PAVSS) . Each  system  consists  of 
sensors  (propeller  anemometers  or  acoustic  radarsj  strategically 
placed  on  lines  perpendicular  to  the  extended  runway  centerline. 

The  instrumentation  (two  GWVSS  lines  and  one  PAVSS  line)  was  set 
up  in  the  middle-marker  to  runway- threshold  region  of  runway  31R. 

In  addition,  a meteorological  tower  was  erected  to  monitor  simul- 
taneously the  winds.  Reference  2 describes  the  test  site  and  the 
operation  of  the  vortex-sensing  systems. 

At  the  time  the  Kennedy  site  was  being  established,  the  Air 
Traffic  Service  of  the  FAA  requested  that  the  Instrument  Flight 
Rules  (IFR)  and  vortex  separations  applied  under  TFR  be  checked  to 
ascertain  their  adequacy  for  protecting  airliners  from  an  inadver- 
tent vortex  encounter.  Because  of  the  expected  limited  traffic 
using  runway  31R,  it  was  decided  to  instrument  a second  runway. 
Since  Kennedy  did  not  have  adequate  real  estate  for  the  vortex- 
tracking equipment  in  the  middle-marker  area  of  the  other  runways, 


a second  airport  was  selected- -Denver ' s Stapleton  International 
Airport.  The  choice  of  Stapleton  was  prompted  by:  (1)  available 
real  estate  in  the  middle-marker  region  of  runway  26L;  ( 2 J the 
different  climate  of  Denver,  Colorado,  as  compared  with  Jamaica, 

New  York;  and  (3)  the  chance  to  record  aircraft  types  which  do  not 
frequent  Kennedy  (e.g.,  the  B-737).  Two  CWVSS,  one  PAVSS,  and  two 
meteorological  towers  were  set  up  at  Stapleton  (Ref.  2). 

Data  collection  began  at  Kennedy  in  July  1973,  and  at  Staple- 
ton  in  August  1973.  By  November  1973,  at  which  time  the  Stapleton 
site  was  closed,  vortices  from  10,000  aircraft  had  been  monitored 
(3100  at  Kennedy  and  6900  at  Stapleton).  The  data  demonstrated 
that  under  IPR,  current  vortex  separations  were  indeed  adequate 
for  preventing  a vortex  encounter.  In  fact,  the  data  indicated 
that  for  most  of  the  time  the  separations  were  unnecessarily  re- 
strictive (Ref.  2). 

For  the  next  2-h  years,  the  Kennedy  site  was  used  as  it  was 
originally  intended--an  operational  site  for  testing  and  evaluat- 
ing vortex  sensors. 

Since  June  1970,  almost  400  unsolicited  reports  of  possible 
vortex  encounters  (incidents)  were  received  by  air  traffic  control 
officers  at  the  London  Heathrow  International  Airport.  Conse- 
quently, the  British  Civil  Aviation  Authority  instituted  a program 
to  gather  information  on  vortex  behavior  under  operational  condi- 
tions, and  on  the  effect  of  a wake  vortex  encounter  by  civil  air- 
craft (Ref.  6).  The  majority  of  reported  incidents  occurred  at  or 
near  Heathrow  and  on  final  approach.  No  accidents  were  recorded 
however.  In  some  cases,  the  encounter  took  place  very  near  to  the 
ground.  Incidents  were  reported  from  a wide  variety  of  aircraft 
pairs  ( leader/fol lower ) , and  it  was  found  that  the  heaviest  jets 
( B - 7 4 7 and  L-1011)  were  cited  in  40  percent  of  all  reported  inci- 
dents, even  though  they  constituted  only  about  12  percent  of  all 
traffic  at  Heathrow  during  peak  periods  (Refs.  7 and  8).  After 
consideration  of  the  incident  reports  and  consultation  with  the 
appropriate  operations  groups,  the  approach-separation  distance  for 
lighter  aircraft  following  a wide-body  jet  in  the  United  Kingdom 


was  increased  in  March  1974  from  five  to  six  nautical  miles. 


There  has  been  close  liaison  between  the  British  Civil  Avia- 
tion Authority  and  the  BAA  on  wake  vortex  research  for  some  years. 
They  jointly  agreed  in  late  1975  that  it  would  be  beneficial  if 
equipment  similar  to  that  tested  at  Stapleton  and  Kennedy  were  in- 
stalled at  Heathrow.  The  test  program  would  afford  the  opportunity 
to  expand  substantially  the  vortex  track  and  meteorological  data 
base  under  new  and  varied  environmental  conditions,  to  correlate 
reported  vortex  incidents  with  measured  vortex  and  meteorological 
conditions,  and  to  track  vortices  from  aircraft  rarely  seen  in  the 
United  States  (e.g.,  Trident,  Viscount,  and  A-300) . 

The  equipment  (two  GIWSS,  one  PAVSS,  and  two  meteorological 
towers)  was  emplaced  between  the  middle-marker  and  runway- threshold 
area  of  runway  28R,  and  became  fully  operational  in  May  1974.  Data 
collection  continued  through  June  1975.  A total  of  12,950  landings 
was  monitored  (Refs.  4 and  5). 

As  a result  of  the  Heathrow  tests,  the  correlation  between  the 
ambient  meteorology  and  vortex  behavior  was  underscored  (Refs.  4 
and  5).  It  was  noted  that  it  would  be  safe  to  use  decreased  sepa- 
rations often;  thus,  the  concept  of  the  VAS  evolved.  The  concept 
of  the  VAS  consists  of  measuring  the  wind  in  the  approach  region, 
comparing  the  wind  velocity  with  a wind  criterion,  and  indicating 
to  the  air  traffic  controller  when  separations  could  be  safely  de- 
creased . 

Tests  were  continued  at  Kennedy  to  assist  in  the  design  of 
the  VAS  and  to  study  the  decay  of  vortices  in  the  terminal  area. 
Portions  of  the  tests  began  in  the  fall  of  1975,  but  the  full  tests 
began  in  March  1976  and  continued  to  the  close  of  the  Kennedy  site 
in  January  1977.  Over  4700  aircraft  passages  were  monitored  in  the 
March  1976  to  January  1977  time  frame. 

The  decay  of  vortices  was  studied  at  Kennedy  using  a Mono- 
static Acoustic  Vortex  - Sens i ng  System  (MAVSS)  (Refs.  9 and  10). 

The  analysis  of  the  data  led  to  a revision  of  the  separation 
standards  for  aircraft  with  a maximum  certificated  gross  takeoff 
weight  of  less  than,  or  equal  to,  12,500  pounds  (classified  as 
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Small  aircraft  by  Air  Traffic  Control)  following  Large  (Large-- 
aircraft  of  more  than  12,500  pounds  maximum  certificated  takeoff 
weight  up  to  300,000  pounds)  and  Heavy  aircraft.  The  new  separa- 
tion rules,  applied  at  the  runway  threshold,  were  promulgated  in 
November  1975;  few  vortex- related  accidents  have  occurred  since 
then.  The  present  ITR  landing  separation  requirements  are:  three 
nautical  miles  for  any  aircraft  following  a Small  aircraft  and  for 
Large  and  Heavy  aircraft  following  a Large  aircraft;  four  nautical 
miles  for  a Small  following  a Large  and  for  a Heavy  following 
another  Heavy;  five  nautical  miles  for  a Large  following  a Heavy; 
and  six  nautical  miles  for  a Small  following  a Heavy. 

The  primary  purpose  of  the  tests  at  the  Chicago  O'Hare  Inter- 
national Airport  was  the  evaluation  of  the  VAS  concept.  The  mid- 
dle-marker region  of  runways  321.,  14R,  and  27R  were  each  instru- 
mented with  a GK'VSS  line.  Outputs  from  the  VAS  were  compared  with 
actual  vortex  behavior  measured  by  the  three  GIVVSS  lines.  Between 
July  1976  and  September  1977,  over  22,500  aircraft  landings  were 
monitored.  In  September  1977,  the  development  tests  of  the  VAS 
were  terminated;  an  upgraded  VAS  for  operational  use  is  now  being 
implemented  at  O'Hare. 

Analysis  of  vortex  behavior  from  over  50,000  aircraft  landings 
showed  that  a wind-rose  criterion  could  be  used  to  determine  when 
separations  may  be  reduced  uniformly  to  three  nautical  miles  be- 
tween the  middle  marker  and  threshold  regardless  of  leader/fol lower 
aircraft  type  (Refs.  2,  4,  and  5).  Most  of  the  time,  vortices  in 
the  middle-marker  to  touch-down  region  either  transport  away  from 
the  extended  runway  centerline  or  decay  to  an  innocuous  level. 

The  few  cases  when  a vortex  persisted  near  the  extended  runway  cen- 
terline for  a time  commensurate  with,  or  in  excess  of  the  three- 
nautical-mile  standard,  were  studied  in  great  detail.  It  was  ob- 
served that  the  measured  one-minute  averaged  winds  could  be  used 
to  predict  when  vortices  would  not  persist  near  the  extended  runway 
centerline.  Whenever  the  wind  exceeds  the  wind-rose  criterion, 
uniform  three-nautical-mile  spacings  may  he  used  with  safety  vor- 
texwise  inside  the  middle  marker  (the  region  where  the  vortex  data 
were  collected).  Not  one  of  the  vortices  from  the  more  than 
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50,000  aircraft  would  have  posed  a aft\  problem  inside  the  middle 
marker  if  all  interarrival  spacing*  w,-.  three  nautical  miles  dur- 
ing the  times  when  the  wind  exceeded  the  wind-rose  criterion. 
Therefore,  it  is  asserted  that  the  safety  of  the  VAS  has  been  dem- 
onstrated for  the  middle-marker  to  runway- threshold  region.  For 
details  about  this  data  collection  and  analyses,  see  Refs.  2,  4, 
and  5. 

3.2  VORTEX  ADVISORY  SYSTEM  DESIGN 

The  VAS  was  designed  to  use  a wind-rose  criterion.  The  system 
compares  the  measured  wind  magnitude  and  direction  (with  respect  to 
each  runway  heading)  with  the  wind  criterion.  The  result  of  the 
comparison  is  indicated  via  a simple  display;  a green  light  means 
that  three-nautical-mile  spacings  may  be  used  for  landing  aircraft; 
and  a red  light  means  that  the  normal  three,  four,  five,  or  six-nau- 
tical-mile spacing  should  be  used  under  IFR  depending  upon  the  re- 
spective aircraft  types.  No  vortex  tracking  sensors  are  used. 

The  VAS  consists  of  four  major  subsystems;  (1)  a meteoro- 
logical subsystem  for  the  measurement  of  the  winds;  (2)  a data- 
processing  subsystem  for  the  processing  of  the  wind  data,  and, 
using  the  VAS  algorithm,  for  determining  when  spacings  between 
interarrival  aircraft  may  be  reduced;  (3)  a data-display  subsystem 
for  the  display  of  separation  requirements  and  wind  conditions  to 
the  air  traffic  controllers;  and  (4)  a performance-monitoring  and 
data- recording  subsystem  for  maintenance  and  archival  purposes. 

3.2.1  Meteorological  Subsystem 

The  meteorological  subsystem  consists  of  a network  of  instru- 
mented towers  placed  about  the  airport  perimeter.  In  concept,  each 
runway  end  would  have  a single  50- foot  tower  approximately  halfway 
between  the  runway  threshold  and  the  middle  marker  and  about  1000 
feet  to  one  side  (to  prevent  vortex  impingement  on  the  tower  dis- 
turbing the  meteorological  measurements).  The  proximity  of  runway 
thresholds,  however,  can  often  permit  the  placement  of  a single 
tower  to  serve  two  (or  more)  runways.  Seven  towers  are  used  in 
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the  O'Hare  system  to  provide  wind  data  on  12  runway  ends. 


fach  50- foot  tower  is  instrumented  witli  3 sets  of  wind  magni- 
tude and  j.1  i rev t ion  sensors,  one  sensor  at  the  50- foot  height  and 
the  remaining  two  at  47  feet.  The  redundancy  provided  by  a triple- 
sensor installation  greatly  increases  system  reliability,  insuring 
acquisition  of  valid  wind  data  and  detection  of  sensor  failures. 

The  instrumentation  transmitting  the  meteorological  data  from 
each  tower  to  a central  facility  consists  of  a multiplexer  which 
sequentially  samples  the  sensor  outputs,  and  a line  modem  which 
serializes  the  data  and  transmits  them  over  a wire  pair  to  receiv- 
ers in  the  control  tower.  A lb-channel,  12-bit  multiplexer  is 
used.  The  multiplexer  operates  under  the  control  of  the  modem 
which  commands  the  scan  rate.  The  modem  operates  in  a lines- 
switching  mode  at  a crystal  controlled  5440-Hz  bit  rate.  In  addi- 
tion to  the  six  multiplexer  channels  used  to  read  the  three  wind 
speed-and-direction  outputs,  four  channels  are  used  to  monitor  the 
status  of  the  tower  electronics  by  monitoring  a precision-voltage 
reference  and  power-supply  outputs. 

Tower  electronics  are  housed  in  an  environmental  enclosure 
mounted  near  the  base  of  each  tower.  Since  lightning  strikes  are 
a major  problem  in  this  type  of  installation,  extra  care  was  taken 
to  insure  against  lightning  damage.  All  input  and  output  signal 
lines  are  protected  with  transient  arrestors.  The  input  60-Hz 
power  line  is  regulated,  and  contains  a separate  transient  arrest- 
or. Standard  FAA  control  lines  are  used  to  transmit  the  data  from 
each  meteorological  tower  to  the  control  tower. 

3.2.2  Pa t a - process ing  Subsystem 

The  serial  data  stream  from  each  meteorological  tower  is  re- 
ceived by  a modem  which  converts  the  input  into  parallel  16-bit 
words  representing  the  output  of  each  channel  sampled  by  the  tower 
instrumentation.  The  output  from  each  receiving  modem  is  input  to 
a microprocessor  (one  for  each  meteorological  tower).  The  micro- 
processors sample  the  wind  data  at  a two  samples/second  rate.  The 
sampled  wind  magnitude  f R ^ ) and  wind  direction  (0 . ) are  used  to 


compute  a one-minute  running  average  (R  and  0)  by  the  following 
scheme:  for  each  sample,  compute  U.  = R.sintj  and  Vj  = R^cosSj; 

then,  compute  U and  V using  a running  128-sample  average,  and  com- 
pute R = (U2  + V 2 j and  0 = tan  ^ (V/U) . A 1 -minute  average  (ac- 
tually a 64-second  average)  was  chosen  since  the  average  life  of  a 
vortex  was  found  to  he  1 minute  (Refs.  2 and  5). 

The  microprocessor  also  performs  the  functions  of  failure 
detection  and  gust  computation.  The  sampled  R^  and  0.  from  each 
sensor  on  a tower  are  compared  at  the  end  of  each  sampling  interval 
(h  second),  and  must  agree  to  within  3 knots  and  20  degrees,  re- 
spectively; if  they  do  not  agree,  at  least  one  of  the  sensors  is 
assumed  to  have  failed.  Normally,  the  50-foot  sensor  data  are 
selected;  if  a 50-foot  R^  or  0.  fails,  the  microprocessor  switches 
to  the  47-foot  Rj  or  0^  which  is  not  in  the  wind  shadow  of  the 
tower.  Failure  of  at  least  two  R^'s  or  B.'s  to  agree  for  eight 
successive  samples  (4  seconds)  causes  a tower- f a i lure  signal  to  be 
generated.  The  microprocessor  calculates  the  wind-gust  magnitude 
using  a sliding  32-second  interval.  Within  each  32-second  interval 
the  sampled  wind  magnitude  is  averaged  using  a 4-sample  running 
average.  Momentary  peaks  due  to  high-frequency  gusts,  which  would 
not  affect  aircraft  operations,  are  filtered  out  by  the  4-sample 
running  average.  Any  measured  peak  must  be  at  least  9 knots  above 
R to  be  considered  a gust,  and  the  gust  value  is  the  peak  value 
observed  during  each  sliding  32-second  interval. 

In  addition  to  outputs  of  R,  9,  and  gust  (if  any),  the  micro- 
processor outputs  system  status  words  to  indicate  which  specific 
failure  (if  any)  is  detected.  Failure  indications  are  displayed 
on  the  system-maintenance  console,  thereby  providing  maintenance 
personnel  with  the  means  to  effect  rapid  repairs. 

The  microprocessor  also  contains  the  VAS  wind-criterion 
algorithm  in  look-up  table  form  for  determining  the  separation 
standard:  the  3/4/5/6  nautical  miles  or  3 nautical  miles  for  all 
aircraft.  The  VAS  wind  criterion  is  the  inner  of  the  two  concen- 
tric ellipses  in  Fig.  5.  The  inner  ellipse  has  major  and  minor 
axes  of  12.5  and  5.5  knots,  respectively;  the  outer  ellipse,  14.5 
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and  7.5  knots,  respectively.  The  major  axes  are  aligned  in  the 
direction  of  the  runway.  If  the  averaged  wind  vector  is  on  or  in- 
side the  inner  ellipse,  the  3/4/5/6-nautical -mile  separations  are 
indicated.  If  the  averaged  wind  vector  is  on  or  outside  the  outer 
ellipse,  the  uniform  3- naut i ca 1 - mi le  separation  is  indicated.  The 
region  between  the  two  ellipses  serves  as  a buffer  zone  to  prevent 
rapid  changes  between  the  3/4/5/ 6-  and  the  3-nautical -mi le  separa- 
tion indications.  If  the  averaged  wind  vector  is  on  or  inside  the 
inner  ellipse,  the  wind  must  increase  so  as  to  reach  the  outer 
ellipse  before  the  3-nautical-mile  separations  are  indicated.  If 
the  averaged  wind  vector  is  on  or  outside  the  outer  ellipse,  the 
wind  must  decrease  so  as  to  reach  the  inner  ellipse  before  the 
3/4/5/6-nautical- mile  separations  are  indicated. 

3.2.3  Data-Display  Subsystem 

Two  types  of  displays  are  used  in  the  VAS,  a system-monitor 
display  and  a runway-monitor  display. 

3. 2. 3.1  System-Monitor  Display 

The  system-monitor  display  is  intended  for  use  by  the  tower- 
| cab  and  IFR-room  supervisors.  The  display  indicates  in  summary 

form  the  winds  measured  by  all  the  meteorological  towers.  The 
primary  function  of  the  display  is  to  provide  an  overview  of  the 
wind  conditions  across  the  airport  enabling  the  supervisor  to 
select  an  operating  configuration  which  will  maximize  traffic  flow. 

3.2.3. 2.  Runway-Monitor  Display 

The  runway-monitor  display  is  intended  for  use  by  a controller 
responsible  for  traffic  on  a single  runway.  The  controller  selects 
the  specific  runway  via  a set  of  thumbwheel  switches.  The  control- 
ler also  indicates  if  arrival  or  departure  winds  are  desired,  e.g., 

I enter  A521.  for  arrivals  on  runway  32L,  and  D32L  for  departures 

from  32b.  The  display  thereafter  accepts  data  with  the  correspond- 
ing label  from  the  data  bus.  Thus,  if  A32L  is  entered,  wind  param- 
eters measured  at  the  tower  near  the  approach  end  of  runway  32L 
are  displayed,  while  a D32L  entry  causes  wind  parameters  measured 
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at  the  tower  near  the  approach  end  of  runway  1 4 R to  he  displayed. 

Separations  are  indicated  by  a red  or  green  light  and  are 
indicated  only  when  an  arrival  runway  is  selected.  A red  light 
means  that  3/4/ 5/6- naut ica 1 - mi le  separations  need  to  be  maintained 
during  IFR;  a green  light  means  that  all  3-nautical-mile  separa- 
tions may  be  used.  The  dimensions  of  the  VAS  wind  criterion 
(inner  ellipse  of  Fig.  5)  have  been  defined  to  permit  aircraft  at 
or  within  the  middle  marker  to  safely  land  if  the  VAS  suddenly 
transitions  from  green  to  red. 

3.2.4  Performance-monitoring  and  Da  t a - record  i ne  Subsystem 


To  facilitate  maintenance  of  the  VAS,  a system  maintenance 
panel  indicates  the  status  of  the  various  components  and  if  re- 
pairs are  needed.  All  the  data  output  by  the  VAS  are  recorded  on 
nine-track  digital  tape.  This  tape  contains  a complete  record  of 
all  VAS  operations  for  use  in  system  diagnostics  and  to  meet  FAA 
operational  requirements  for  a record  of  all  air  traffic  control 
ope  rations. 


4.  VORTEX  ENCOUNTER  MODEL 


The  wind  criterion  employed  by  the  VAS  ensures  that  vortices 
will  not  pose  a problem  to  a landing  aircraft  between  the  middle 
marker  (MM)  and  touch  down  when  the  green  light  is  on.  Permitting 
separations  to  decrease  below  the  4-,  5-,  and  6-naut ical -mi le 
standards  inside  the  MM  will  contribute  to  minimizing  delays.  How- 
ever, significant  improvements  in  minimizing  delays  can  be  realized 
if  the  coverage  of  the  VAS-protected  region  can  be  extended  to  the 
outer  marker  (OM) . Measurements  have  been  made  of  vortex  behavior 
between  the  MM  and  OM;  Volume  II  of  this  report  addresses  a com- 
prehensive vortex  data  collection  program  in  this  region. 

The  objectives  of  this  section  are  to  derive  a general  model 
for  calculating  the  probability  of  a hazardous  vortex  encounter 
between  the  OM  and  the  MM  and  to  determine  the  applicability  of  the 
VAS  green  light  to  this  region.  The  intent  is  to  be  as  general  as 
possible,  so  that  the  results  may  be  applied  at  any  airport. 
Limitations  of  the  model  (i.e.,  anomalous  vortex  behavior,  terrain 
effects,  etc.)  will  be  delineated  herein  or  in  the  appendices. 
However,  the  model  is  based  on  state-of-the-art  knowledge  of  how 
vortices,  aircraft,  and  winds  behave. 

Two  independent  conditions  must  be  met  for  a hazardous  vortex 
encounter  to  occur.  First,  the  trailing  vortex  system  must  be 
sufficiently  strong  to  be  hazardous  to  the  vortex-encountering  or 
following  aircraft,  and  second,  the  following  aircraft  must  encounter 
the  vortices  of  the  leading  or  vortex-generating  aircraft.  Thus, 
two  probabilities  are  defined,  the  probability  that  a vortex  is 
hazardous,  P^,  which  is  related  to  the  strength  of  the  trailing 
vortex  system  and  to  the  size  and  control  characteristics  of  the 
vortex -encounter ing  aircraft;  and  the  probability  of  encountering 
the  vortices,  P , which  is  related  to  the  flight  paths  of  the  two 
aircraft  and  to  the  motion  of  the  vortex  system.  The  probability 
of  a hazardous  vortex  encounter,  Pj  , is  then 


4.1  PROBABILITY  OF  A VORTliX  P.NCOUNTLR 

4.1.1  Oef_i  fi  i t i oi!  of  Vortex  Encounter 

Before  deriving  an  analytic  expression  for  the  probability  of 
a hazardous  vortex  encounter,  it  is  necessary  to  define  what  is 
meant  by  an  encounter.  Although  the  definition  is  simple  and  by 
no  means  complete  or  unique,  it  docs  take  into  account  a practical 
criterion  for  a hazardous  vortex  encounter  to  occur.  This  cri- 
terion is  that  the  swirling  velocities  of  the  vortex  induce  veloc- 
ities on  the  wing  of  the  vortex-encountering  aircraft  in  excess  of 
the  roll  velocities  which  result  with  full  aileron  deflection.  In 
other  words,  a hazardous  condition  exists  when  the  vortex  induces 
roll  rates  exceeding  the  roll  control  authority  of  the  vortex-pene- 
trating aircraft.  Analysis  of  flight  test  (Refs.  11  and  12)  and 
simulation  (Refs.  13  and  14)  results  has  verified  that  the  cri- 
terion agrees  with  pilot  assessment  of  a hazardous  vortex  en- 
counter. To  be  conserrat i ve  in  the  model,  a more  stringent  cri- 
terion shall  be  used:  a hazardous  condition  is  defined  to  exist 
when  the  vortex- induced  roll  rate  exceeds  some  fraction  f (to  be 
defined,  but  less  than  1)  of  the  roll  control  authority  of  the 
vortex- encounter ing  aircraft.  Since  the  aircraft  are  in  trail, 
only  the  worst  case  of  an  axial  vortex  encounter  is  considered 
herein . 

The  vortex  wake  at  a distance  behind  an  aircraft  is,  in 
general,  comprised  of  two  counterrotating  vortices  spaced  a dis- 
tance bQ  apart.  Sufficiently  far  from  the  center  of  each  vortex, 
the  swirling  velocity  v is  described  by 


where  T is  the  circulation  or  strength  of  the  vortex,  and  R is  the 
distance  from  the  center  of  the  vortex.  Inside  the  circle  of 
radius  R,  the  swirling  velocities  are  larger  than  v although  the 
swirling  velocity  is  zero  at  the  center  of  the  vortex.  Define  Rq 
as  the  radius  at  which  the  swirl  velocity  equals  the  fraction  f of 


the  maximum  ai leron- induced  wingtip  velocity  v.j,: 


R 
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Rq  defines  a cross  section,  such  that  if  the  wingspan  of  an  air- 
craft is  embedded  in  the  circle  of  radius  Rq  , the  fraction  f of 
the  roll-control  authority  of  the  aircraft  is  exceeded.  Thus,  for 
a potentially  hazardous  vortex  encounter  the  vortex  tangential 
velocity  at  Ro  must  be  greater  than  or  equal  to  fv,j.. 

For  a vortex- encounter ing  aircraft  of  maximum  roll  capability 
p and  wingspan  b , 


where  pbc/2U  = p is  the  maximum  non-dimensional  roll  rate  and  U is 
the  airspeed  of  the  vortex-encountering  aircraft.  Thus, 


0 2 m fpU 

defines  a hazard  radius  about  each  vortex  dependent  on  vortex 
strength,  maximum  non-dimensional  roll  rate,  fraction  of  the  maxi- 
mum ai leron- induced  wingtip  velocity  imparted  to  counter  the  vor- 
tex-swirling velocity,  and  airspeed. 

For  commercial  jetliners,  the  maximum  tip  velocity  for  full 
aileron  deflection  is, at  a minimum,  six  percent  of  the  flight 
speed;  for  general-aviation- type  aircraft,  the  maximum  tip  velocity 
is  at  least  eight  percent  of  the  flight  speed.  (The  maximum  non- 
dimensional  roll  rate  is  the  product  of  the  a i leron- and- spo i ler 
control  effectiveness  and  the  maximum  a i leron- and-spoi ler  deflec- 
tion angles  divided  by  the  roll-damping  coefficient  (Ref.  15). 
Ailerons  and  spoilers  for  new  commercial  jetliners  are  designed  so 
that  p is  at  least  0.06.) 

Inside  a circle  of  radius  R , the  characteristic  vortex- 

o 


swirling  velocities  are  larger  than  the  maximum  wingtip  velocities 
attainable  with  the  fraction  f of  the  roll  capability  of  the  vor- 
tex-encountering aircraft.  If  the  wingspan  of  the  vortex-encoun- 
tering aircraft  is  smaller  than  2Rq,  the  aircraft  can  be  entirely 
immersed  in  a swirling  flow  whose  tangential  velocity  exceeds  the 
fraction  f of  the  aircraft's  roll  authority.  On  the  other  hand, 
if  the  wingspan  is  larger  than  2Rq , vortex  - induced  swirl  velocities 
cannot  exceed  the  fraction  f of  the  roll  authority  of  the  aircraft. 
Thus,  the  condition  for  a potentially  hazardous  vortex  encounter 
is  that 


R > b /2. 
o e 


(6) 


According  to  the  model  then,  a potentially  hazardous  vortex 


encounter  can  occur  only  if  R 


> be/2, 


and  if  a following  aircraft 


is  fully  located  within  the  circle  defined  by  Rq  . The  probability 
of  a potential  hazard  is: 


( 1 R > b /2, 

1 o e 

( 0 R < b /2. 

' o e 


(7) 


The  encounter  is  either  hazardous  (P^  = 1)  or  it  is  not  (P^  = 0). 

To  be  conservative  "soft"  or  partial  encounters  are  considered  to 

be  hazardous.  = 0 includes  encounters  of  little  or  no  concern, 

and  P^  = 1 includes  cases  which  are  serious  as  well  as  many  which 

are  of  little  concern.  Thus,  P,  = P,P  is  defined  as  zero  when 

’he  he 

the  hazard  radius  Rq  decays  below  be/2,  and  is  finite  for  a spe- 
cific following  aircraft  when  Rq  exceeds  be/2.  Note  that  P^e  = 0 
does  not  preclude  a vortex  encounter,  it  means  that  an  encounter 
is  not  hazardous. 

4.1.2  Vortex  Decay 

The  hazard  posed  by  a vortex  primarily  involves  the  vortex 
strength.  The  strength  or  circulation  of  a vortex  decays  with 
distance  behind  the  vortex  - gencrat i ng  aircraft,  and  the  decay  is 
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influenced  by  many  factors.  Close  behind  the  aircraft,  the  ini- 
tial dispersion  of  the  wake  is  governed  by  configuration,  which 
includes  geometry  (particularly  flap  and  engine  placement)  and 
thrust  and  flap  settings. 

While  the  initial  vortex-decay  rate  is  configuration-sensi- 
tive, ultimately  the  decay  rate  must  he  determined  by  self-decay 
of  the  vortex  pair  coupled  with  atmospheric  shear,  turbulence,  and 
stratification.  While  there  is  some  understanding  of  the  turbulent 
mechanisms  involved  in  dissipation  of  vortices,  no  theories  which 
can  reliably  predict  the  decay  of  an  aircraft  wake  are  available 
today.  Theories  do  exist  for  describing  the  onset  of  either  of 
the  two  catastrophic  decay  modes:  sinuous  instability  or  vortex- 
linking (Ref.  16)  and  vortex-bursting  (Ref.  17).  The  predicted 
times  are  generally  within  a factor  of  two  of  the  observed  times 
for  both  breakup  hy  linking  and  by  bursting.  Atmospheric  turbu 
lence  is  the  mechanism  that  causes  catastrophic  decay  to  occur. 

To  be  conservative  in  the  model,  the  catastrophic  decay  modes 
(which  occur  most  of  the  time)  will  be  ignored;  only  the  slowest 
decay  mode  (viscous  decay)  will  be  considered. 

Decay  by  eddy  viscosity  involves  extremely  calm  atmospheric 
conditions,  when  linking  or  bursting  may  not  be  occurring,  but 
there  is  still  sufficient  aircraft-generated  turbulence  to  cause 
eventual  vortex  erosion.  Aircraft  in  a "dirty"  configuration 
(gear  and  flaps  down)  generate  turbulence  which  is  added  to  the 
wake  permitting  eddy-viscosity  decay  of  the  wake.  Considering 
only  this  mode  of  decay  establishes  an  outer  limit  or  upper  bound 
on  the  life  of  a vortex. 

Eddy- viscosity-caused  decay  of  the  vortices  behind  "dirty" 
aircraft  is  illustrated  in  Figure  6 (Ref.  18).  The  sources  of  the 
data  are  the  ground-based  measurements  taken  by  TSC  at  Kennedy 
International  Airport  (Ref.  10),  NASA/FAA  flight  test  measurements 
using  probing  aircraft  (Fig.  25  of  Ref.  19),  and  NASA  towing-tank 
results  (Ref.  20).  Although  the  definition  of  vortex  strength 
differed  for  the  different  data  sources,  the  major  results  would 
not  differ  if  other  definitions  were  used. 


The  figure  shows  that  non- d imens i ona 1 i z i ng  the  circulation 
by  an  effective  initial  circulation  fQ  and  non-dimensional izing 
distance  d behind  the  vo r t ex  - gene ra t i ng  aircraft  by  its  wingspan 


b , i t : 
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ift  coefficient  , and  its  aspect  ratio  A,  cause  the 


data  points  to  coalesce  fairly  well.  It  is  remarkable  that  such 
diverse  measurements  indicate  similar  decay.  Beyond  d(^/b  A of 
9.58  (corresponding  to  the  region  beyond  about  50  span  lengths), 
the  data  fit  is 


~ = 9.58 
o 


(W1 


(8) 


In  this  region,  the  eddy  viscosity  acts  to  diffuse  the  vorticity 
causing  wake  dissipation  to  proceed  as  d or  (time)  . Earlier, 
when  r/r0  is  essentially  constant,  aircraft  configuration  effects 
are  acting  to  diffuse  vorticity  so  that  decay  can  subsequently 
occur . 

The  shape  of  the  curve  in  Figure  6 showing  constant  F at  short 
distances  and  thereafter  a (time)  1 dependence  is  the  same  as  found 
for  the  core  circulation  of  a B-747  in  landing  configuration 
measured  by  a laser  Doppler  velocimeter  (Fig.  5 in  Ref.  21).  The 
"break"  for  the  laser  measurements  occurred  at  about  33  span 
lengths,  somewhat  closer  to  the  vortex-generating  aircraft  than 
the  "break"  shown  in  Figure  6.  Other  investigators  (several 
papers  in  Ref.  22)  have  found  slightly  different  shaped  curves  or 
noted  decay  to  be  stronger  near  the  ground,  but  the  main  feature 
of  constant  circulation  followed  by  decay  appears  consistent  with 
al 1 observat ions . 

Figure  6 shows  that  circulation  is  reduced  to  20  percent  of 
its  initial  magnitude  by  dCj /b  A = 50,  about  250  span  lengths  back 
( A / C | is  about  5 for  most  landing  aircraft  during  the  approach). 

For  a B-747  at  approach  speed,  this  is  about  8 nautical  miles  back 
at  a time  of  205  seconds.  For  a B-727  at  approach  speed,  this  is 
about  4.4  nautical  miles  back  at  a time  of  131  seconds.  If  the 
ambient  turbulence  is  light  or  greater,  catastrophic  decay  through 
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linking  or  core-bursting  will  occur  before  these  times.  However, 
only  eddy-viscosity  decay  will  be  assumed  in  the  model  to  insure 
that  the  model  is  conservative.  The  hazard  radius,  R , thus 
becomes 


R = 
o 
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As  the  distance  behind  the  vor tex - generat i ng  aircraft,  d,  increases, 
the  hazard  radius  Rq  decreases.  Recalling  equation  (6J, 


7.6  T b 
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is  required  for  a potentially  hazardous  vortex  encounter.  Thus, 
= 1 fand  Pj  > 0)  only  when 


15.2  r. 
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for  9.58. 
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(11) 
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When  the  solution  to  equation  (11)  yields  d <9.58  (C./b  A) 

*-*  s 

then  R can  never  exceed  b_/2  at  any  distance,  thus  there  is  no 
o & 

distance  at  which  a potentially  hazardous  encounter  can  occur,  and 


d = 0. 


(12) 


Measurements  (Refs.  2,  21,  and  various  papers  in  22)  show 
that  the  circulation  or  strength  of  a vortex  increases  with  the 
spatial  extent  about  a vortex  up  to  an  ultimate  value  ; the 
radius  at  which  r is  attained  is  denoted  by  r . If  the  vortex 
field  is  sampled  (such  as  by  an  aircraft  wing)  with  a probe  with 
r<roo  , only  a portion  of  will  be  experienced.  Given  two  air- 
craft with  different  wingspans,  as  long  as  b^  < 2roj,  the  strength 
of  the  vortex  experienced  by  the  aircraft  with  the  larger  wingspan 
will  be  greater  than  for  the  smaller  wingspan.  Various  radial 
dependencies  have  been  claimed  by  researchers  (Ref.  2);  a linear 
variation  with  be  for  r is  adopted  herein  as  it  is  consistent  with 
strength  data  measured  using  real  aircraft  (Refs.  10,  21,  and  23) 
rather  than  aircraft  models  in  towing  tanks  (Refs.  2 and  17). 

liquation  (11)  requires  an  estimate  of  the  initial  strength 

T of  a vortex  (strength  before  the  onset  of  decay).  The  initial 
o 

strength  Fo  of  the  vortices  as  would  be  experienced  by  an  encoun- 
tering aircraft  of  wingspan  bp  is  shown  in  Appendix  A to  have  the 
form 


where  m and  r.  are  two  curve- f it  ting  parameters;  values  of  m and 
Fjnt  are  tabulated  in  Appendix  A.  Equation  (11)  may  now  be  re- 
written (dCj/b^A  2l  9.58), 
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15.2 

fpU 


+ r . j . 
int  ’ 


(14) 


When  dGj /h < 9.58,  equation  (12)  is  valid  (i.e.,  no  hazard). 

The  NTSB  data  base  (Section  2)  indicated  that  to  the  extent 
separation  data  were  available  all  the  probable  vortex-related 
landing  accidents  occurred  when  the  separations  were  less  than  the 
respective  separation  standards.  To  determine  an  upper  bound  for 
the  fraction  f of  the  maximum  roll-control  authority,  the  distance 
d was  set  equal  to  the  respective  separation  standard.  The 
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standards  that  pertain  to  the  OM  location  were  used  (3,  4,  and  S 
nautical  miles),  not  those  lor  the  runway  threshold  location  (3,  4, 

5,  and  f>  nautical  miles).  Table  2 shows  the  12  aircraft  con- 
sidered along  with  their  landing  speeds  and  wingspans.  Table  3 
shows  the  upper  bound  on  the  fraction  of  roll-control  authority 
required  to  neutralize  the  vortex-swirl  velocity.  The  largest  f 
is  0.75b,  and  occurred  for  the  PA-28  following  the  DC- 8. 

In  the  spirit  of  being  conservative,  f will  be  set  equal  to 
one -half  the  maximum  value  found  in  Table  3 (f  = 0.5  x (0.756)  = 
0.378).  No  further  justification  is  offered;  an  f = 0.378  means 
that  only  one-half  of  the  maximum  roll  authority  required  to 
safely  fly  the  current  separation  standards  is  expected  to  be 
conservative . 

Using  f = 0.378  and  equation  (14)  (or  equation  (12),  if 
necessary) , the  distances  at  which  the  hazard  radius  becomes 
smaller  than  the  semispan  can  be  calculated.  These  distances  are 
given  in  Table  4.  Ignoring  for  the  moment  the  Sma 1 1 - category  air- 
craft (PA-28  and  Learjet),  the  only  calculated  distances  in  excess 
of  3 nautical  miles  are  for  the  B-727,  B-737,  and  DC-9  following 
the  B-747.  Except  for  these  cases,  separation  standards  could  be 
set  at  3 nautical  miles  between  the  OM  and  MM  with  or  without  a 
VAS . (Recall  that  the  accidents  are  confined  to  the  ground - ef fee t 
region  inside  the  MM.)  Note  that  the  DC-9  must  be  2.54  nautical 
miles  behind  the  DC-10  to  reduce  the  probability  of  hazard  to  zero. 
The  DC-9  accident  at  Dallas-Fort  Worth  Greater  Southwest  Interna- 
tional Airport  occurred  with  the  DC-9  about  2 nautical  miles  be- 
hind the  DC- 1 0 . 

So  far  the  model  has  only  considered  vortex  decay.  To  show 
that  3 nautical  miles  is  safe  behind  the  B-747  when  VAS  is  used, 
the  model  must  now  be  expanded  to  consider  both  vortex  motion  with 
the  wind  and  vortex  descent.  It  is  very  important  to  note,  how- 
ever, that  vortex  decay  alone  does  permit  decreased  separation 
standards  for  most  aircraft  pairs  outside  of  ground  effect. 
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TABLE  2.  AIRCRAFT  PARAMETERS 


Aircraft  Model 

Landing  Speed 

Wingspan 

( ft/  sec.) 

(?t) 

B-707/120B 

232.3 

130.9 

B-727/100 

205 . 8 

108.0 

B-737/100 

197.0 

93.0 

B-  747/200B 

238.0 

195.7 

DC-8/20 

222 . 0 

142.3 

DC-9/20 

189.6 

93.3 

DC-10/30 

232.3 

165.3 

L- 1011/200 

241.1 

155.3 

B-707/32CB 

232.3 

145.8 

DC-8/62 

210.2 

148.4 

PA-28/180 

110.0 

30.0 

Lear j et - 2 5 

154.0 

35.6 

! 


4.1.5  Encounter 


l'robab  i 1 i ty 

The  probability  of  a hazardous  encounter  is  the  product  of  the 
probability  that  the  vortex  is  hazardous,  1’^  (1  i f RQ  ^ be/2),  and 
the  probability  of  encountering  that  hazardous  vortex,  P . This 
section  computes  the  probability  of  encounter  for  the  following 
aircraft  where  encounter  is  defined  as  being  fully  immersed  within 
the  hazard  region  defined  by  RQ . 

figure  7 shows  the  encounter  region  in  the  cross  plane  of  the 
approach  path.  The  origin  is  taken  as  the  nominal  location  of  the 
generator  aircraft  at  the  time  of  its  passage  through  the  cross 
plane.  Normally  this  would  be  the  point  at  which  the  glide  path 
intercepts  the  cross  plane,  as  shown  in  the  figure.  The  lateral 
displacement  from  the  origin  is  denoted  by  y and  vertical  by  z,  as 
shown.  Navigation  errors  of  the  generator  aircraft  are  assumed  to 
be  independent  Normally  distributed  with  standard  deviations 
a and  a ^ , respectively,  statistically  defining  the  generation 
point  of  the  vortex  pair. 

The  follower  aircraft,  as  shown  in  Figure  7 intercepting  the 
cross  plane,  may  be  on  a different  nominal  flight  path  from  that  of 
the  generator.  The  mean  lateral  and  vertical  displacements  of  the 
follower  describing  this  flight  path  are  labeled  L and  G,  re- 
spectively. These  are  also  assumed  independent  Normally  dis- 
tributed with  standard  deviations  a and  a 

ye  z e 

The  position  of  the  center  of  the  vortex  pair  is  a function  of 
descent  rate,  cross  wind,  and  the  time  since  the  passage  of  the 
generator.  IV  and  L)  are  the  resultant  mean  lateral  and  vertical 
transport  components,  with  standard  deviations  and  . (The 
evaluation  of  these  parameters  will  be  considered  further  in  the 
next  sections.)  In  addition  to  the  variance  associated  with  the 
transport  mechanism,  the  position  of  the  vortex  pair  is  also  sub- 
ject to  the  variance  in  the  navigation  of  the  generator,  thus  the 
total  lateral  and  vertical  variances  are: 


(16) 
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FIGURE  7.  ENCOUNTER  GEOMETRY  IN  PLANE 
PERPENDICULAR  TO  FLIGHT  PATH 
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(17) 
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With  the  locations  of  the  center  of  the  vortex  pair  and  the 
follower  aircraft  thus  defined,  the  encounter  can  be  defined  sta- 
tistically and  its  probability  calculated.  Section  4.1.1  defined 
a possibly  hazardous  encounter  to  be  that  condition  where  the  en- 
tire wingspan  of  the  encounter  aircraft  is  within  the  vortex 
hazard  cross  section.  Two  distinct  cases  must  be  evaluated  depend- 
ing upon  the  relationship  of  Rq  to  bQ , the  spacing  between  the  two 

vortices  (b  = n b / 4 ) : 

o g' 

1.  Rq  >>  bQ  --  hazard  radius  much  larger  than  the  separation 
between  the  two  vortices.  In  this  case  the  hazard  cross  sec- 
tion can  be  considered  to  be  a single  circle  of  radius  R 

h o 

located  at  the  center  of  the  vortex  pair. 

2.  R ~ b --  hazard  radius  of  the  same  magnitude  or  smaller 

oo  h 

than  the  separation  between  the  two  vortices.  In  this  case 
the  hazard  cross  section  can  best  be  considered  to  be  two  dis- 
tinct circles  of  smaller  hazard  radius  (corresponding  to  the 
strength  and  circulation  of  the  single  vortex). 

These  two  cases  are  illustrated  in  Figure  8. 

For  the  hazard  potential  to  exist,  the  wingspan  of  the  fol- 
lower must  lie  fully  within  the  R circle.  The  hazardous  encounter 
cross  section  is  defined  by  the  geometric  pattern  mapped  out  by 
fitting  a wingspan  b£  inside  the  circle.  If  2Rq  < b , the  cross 
section  is  zero.  Otherwise,  if  the  center  of  the  encounter  cross 
section  is  located  at  (y  , z ) , the  center  of  the  vortex-cncountcr- 
i ng  aircraft  must  lie  on  the  y-axis  between  y + (Rq  - b e / 2 ) and 
yQ  - (R0  - b e / 2 ) and  on  the  z-axis  between  z^  + ( Rq 2 - b 2/4)1/-  for 
the  aircraft  to  be  completely  engulfed  within  a circle  of  radius 
Rq.  Figure  9 shows  the  geometry.  For  mathematical  convenience, 
let  the  hazardous  encounter  cross  section  be  the  rectangle  (see 
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Figure  9)  of  height  2 f Ro 2 - b 2/4)**  and  width  2(RQ  - be/2)  rather 
than  the  smaller  ellipse-like  shape  (a  conservative  approximat ion). 

The  net  mean  position  of  the  hazardous  encounter  cross  section 
relative  to  the  follower  aircraft  is 

horizontal:  VV  - L 


vertical : D - G . 


The  net  variance  in  that 

2 

horizontal:  Ojl  - 

2 

vertical : Oy  - 


position  is 


2 2 2 

a +o  +o 

yg  ye  w 


2 2 2 

a +o  + on 
z g z e D 


(18) 

(19) 


Thus  the  horizontal  and  vertical  encounter  probabilities  are  the 
probabilities  that  Normal  variates  for  the  aircraft/vortex  spacing 
with  the  above  means  and  variances  lie  within  + (R  - b e/ 2 ) hori- 
zontally and  _+  (R  2 - b e 2 / 4 ) 2 of  zero: 


For  case  (1),  R >>  b , the  probability  that  the  vortex- 
o o 

encountering  aircraft  is  within  the  hazard  rectangle  is  given 
directly  by  the  product  of  P|lf;  and  Py^: 


Pe  ‘ PHE  PVF: 


(22) 
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For  case  (2),  Rq  ~ h , two  smaller  circles  of  radius  RQ//J  are  used 
to  maintain  equal  area  of  the  circles  (and,  hence,  continuity  of 
the  probability)  and  to  reflect  the  effect  and  identity  of  both 
vortices  (Ref.  24).  Assuming  that  both  circles  remain  relatively 
close  to  the  center  of  the  vortex  pair,  then 


HE  VE 


with  P,..,  and  P,.r.  evaluated  at  the  reduced  radius.  The  size  of  R„ 

HE  V E o 

where  case  (1)  transitions  to  case  (2)  is  not  unique;  it  is  merely 
a mathematical  technique  to  avoid  double  counting  areas  of  over- 
lapping hazard  radii  when  R is  small  enough  to  distinguish  the  two 
vortices.  The  break  point  assumed  between  case  (1)  and  case  (2) 
in  the  program  given  in  Appendix  D is  at  Rq  50  percent  above  the 
vortex  separation 


R„  = 1.5  bo 
o o 


The  probability  of  a hazardous  vortex  encounter,  P^  , is  cal- 
culated by  substituting  the  Ro  calculated  [from  equation  (9)]  for 
a given  interaircraft  spacing,  d,  into  the  above  expressions  for 
P.,p  and  P. , and  integrating  over  all  possible  wind  runs  W: 


= P.P 
h e 


p(W)  dW, 


To  determine  the  expected  encounter  risk  for  a given  approach,  Pj 
should  be  integrated  along  the  direction  of  flight  from  the  OM  to 
the  MM.  Instead,  the  maximum  value  of  P^e  for  two  aircraft  at 
minimum  prescribed  separations  is  evaluated  to  determine  the 
maximum  potential  hazard  at  the  minimum  spacing.  The  integration 
along  the  final  is  not  meaningful  as  the  two  aircraft  will  not 
maintain  identical  speeds  or  spacings  between  the  OM  and  MM.  If 
the  following  aircraft  is  slower  (faster)  than  the  leading  aircraft, 


the  separation  will  increase  (decrease)  as  the  aircraft  proceed 

along  the  ILS.  Herein,  the  separation  is  held  fixed,  and  I’^c  is 

evaluated  at  the  location  between  the  OM  and  MM  where  P,  would  be 

he 

at  a maximum  at  the  minimum  legal  spacing  (for  all  aircraft  pairs 
Pjie  is  maximized  at  the  OM  location- -all  calculations  were  done  with 
the  following  aircraft  at  the  OM , and  the  OM  was  assumed  to  be  7 
nautical  miles  from  the  runway  as  this  also  maximizes  P^e) . This 
places  a conservative  upper  bound  on  the  encounter  risk  that  might 
be  experienced  on  final  approach  between  the  outer  and  middle 
ma  rkers . 

4.2  NAVIGATION  MODEL 

It  is  assumed  that  the  navigation  deviations  of  the  vortex- 
generating and  vortex- encounter i ng  aircraft  were  Normally  distrib- 
uted. It  is  further  assumed  that  the  peaks  of  the  lateral  and 
vertical  distributions  are  coincident;  that  is,  the  midpoints  of 
the  two  distributions  intersect  on  the  nominal  flight  path. 

In  a study  conducted  for  the  FAA  (Ref.  25) , an  Unusual  Events 
Recording  System  was  installed  on  a B-737  for  6 months  of  opera- 
tion. Figures  10  and  11  show  the  measured  standard  deviations  for 
IFR  operations.  (The  localizer  course  was  defined  for  a 10,000- 
foot  runway).  3evond  5 nautical  miles,  the  B-737  data  indicate 
extreme  variations  as  the  data  frequently  included  portions  of  a 
turn  preceding  the  localizer  intercept  near  the  OM.  Inside  the 
OM,  the  lateral  deviations  for  the  B-737  IFR  approaches  to  a number 
of  airports  are  the  same  as  for  simultaneous  approaches  to  parallel 
runways  at  Chicago's  O'ilare  International  Airport  by  several  air- 
craft types  (Ref.  26).  Operations  in  the  vertical  plane  covered 
almost  the  entire  + 0. 7-degree  width  of  the  glide-slope  beam  at  a 
range  of  1 to  8 nautical  miles  from  the  runway  threshold.  This 
dispersion  for  IFR  operations  was  greater  than  expected,  but  was 
attributed  to  the  pilot  using  a visual  approach  as  soon  as  the 
airport  was  in  sight. 

The  broken  lines  in  Figures  10  and  11  are  a least  - squares  fit 
to  the  B-737  data  collected  inside  the  OM  (about  26,000  feet  (4.3 
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nautical  miles)  from  the  runway  threshold  for  the  measurements). 
The  equations  for  the  standard  deviation  are: 


o = 0.01121^  - 9.4206, 

y t 

a,  = 0.0039L.  + 9.8049, 


(26) 


where  L.j,  is  the  distance  (in  feet)  from  the  runway  threshold.  The 
correlation  coefficients  of  the  least  - squares  fit  to  a and  o _ are 
0.98  and  0.96,  respectively.  It  is  assumed  that  equations  (26) 
may  be  used  for  any  ITS  runway  between  the  runway  threshold  and 
the  OM  (4  to  7 nautical  miles  from  the  runway  threshold). 


4.5  VORTEX- DESCENT  RATE  MODEL 

To  calculate  the  vertical-encounter  probability  (Eq.  (21)), 
the  vortex-descent  distance  D and  the  standard  deviation  a ^ are 
required.  The  state  of  knowledge  on  vortex  - descent  is  oftentimes 
misinterpreted;  Appendix  B reviews  the  existing  data. 

Vortex-descent  rates  used  herein  were  extracted  from  laser- 
Doppler- velocimeter  measurements  taken  at  the  John  F.  Kennedy 
International  Airport  (Refs.  27-29).  Initial  descent  rates  were 
extracted  from  the  vortex-track  data  given  in  Reference  28;  only 
initial  descent  rates  were  extracted  as  the  vortices  entered  and 
became  dominated  by  ground  effect  within  about  20  seconds.  As  dis- 
cussed in  Reference  30  and  in  Appendix  B,  the  initial  descent  rate 
is  a good  description  of  the  vortex  descent  for  the  first  30  to  60 
seconds  of  the  life  of  a vortex  outside  of  ground  effect.  Cata- 
strophic decay  (leading  to  the  early  demise  of  a vortex)  a. id  at- 
mospheric stability  affect  the  descent  rate,  but  are  neglected. 

Figure  12  shows  B-"47  initial  descent  rates  extracted  from 
the  Kennedy  data.  The  mean  initial  descent  rate  was  found  to  be 
6.3  ft. /sec  with  a standard  deviation  of  1.9  ft/sec.  (For  the 
B-707,  the  mean  initial  descent  rate  was  5.2  ft/scc  with  a standard 
deviation  of  1.8  ft/scc.)  Although  the  standard  deviations  arc 
relatively  large,  most  of  contribution  to  the  large  deviations  came 


FIGURE  12.  INITIAL  DESCENT  RATE  OF  B-747  VORTICES 


from  the  relatively  large  number  of  cases  with  high  descent  rates. 


It  was  not  possible  to  justify  extracting  initial  descent 
rates  for  other  aircraft  types  using  the  Kennedy  data  as  there 
were  too  few  cases  (there  were  at  least  100  cases  each  of  the  B-707 
and  B-747).  Assuming  elliptic  wing  loading,  the  initial  vortex 
circulation  is  given  by  (Ref.  2): 

4 W 

Fo  ~ TrpVb  ’ ( 2 7 1 

g 


where  W is  the  gross  weight  of  the  aircraft,  and  p is  the  density 
of  the  air.  The  theoretical  initial  descent  speed  is  then  given 
by: 


V 


D 


8 W 

7T  pVb  ^ 
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(28) 


(bQ  = ^ b , with  elliptic  loading).  With  p=2. 34x10  3 slug/ff5,  the 
calculated  initial  descent  rates  are  5.3  and  6.8  ft/ sec  for  the 
B-707  and  B-747,  respectively.  (The  values  of  the  parameters  used 
are  given  in  Tables  2 and  5.  The  weights  are  the  maximum  landing 
weights.)  The  calculated  and  measured  initial  descent  rates  agree 
quite  well:  5.3  and  5.2  ft/s ec  for  the  B-707;  6.8  and  6.3  ft/sec 
for  the  B-747.  Thus,  initial  descent  rates  for  other  aircraft 
types  were  calculated  using  the  assumption  of  elliptic  wing  load- 
ing, and  standard  deviations  were  selected  based  on  the  magnitude 
of  the  initial  descent  rates  (see  Table  5). 


4.4  CROSS-WIND  MODEL 

Winds  are  three-dimensional  motions  of  the  air,  and  consist 
of  very  large  to  very  small  scale  temporal  and  spatial  variation  . 
The  variability  of  the  wind  is  caused  and  governed  by  the  rotation 
of  the  earth  (Coriolis  force),  geographic  character! st ics  (such  as 
orographic  effects),  and  the  available  solar  energy  reaching  the 


ieasured  values. 


earth's  atmosphere  and  surface  which  is  further  a function  of  the 
time  of  day.  Other  dominating  factors  causing  wind  variability 
are  land-sea  influences,  terrain  type,  elevation,  available  water, 
vegetation,  and  many  additional  natural  and  artificial  constituents. 

Constructing  a wind  model,  especially  one  which  implies  a 
correlation  between  surface-wind  measurements  and  the  winds  at 
altitudes  above  ground  level  in  excess  of  1000  feet,  is  difficult. 
Adequate  wind  models  exist  for  the  planetary  boundary  layer  (be- 
tween the  ground  and  about  500  feet  in  altitude)  as  the  turning 
of  the  wind  with  height  is  relatively  small.  Above  the  planetary 
boundary  layer,  little  is  known  about  how  to  relate  the  winds  with 
surface  measurements.  However,  the  analysis  of  the  VAS  does  not 
require  a complete  wind  model,  only  one  which  probabilistically 
relates  winds  at  the  surface  in  "green  light"  conditions  to  the 
winds  aloft . 

Cross  winds  play  a role  in  preventing  vortex  encounters  during 
approaches  to  landing.  The  larger  the  cross-wind  component,  the 
less  likely  that  a vortex  encounter  can  occur.  For  small  or  zero 
cross-wind  components,  vortex  encounters  are  prevented  by  the 
dissipation  of  the  vortices,  by  the  vertical  displacement  of  the 
vortices,  and  by  the  relative  flight  paths  of  the  aircraft. 

The  analytical  model  discussed  in  Section  4.1  requires  a 
cross-wind  p.d.f.  to  calculate  the  probability  of  a hazardous- 
vortex  encounter.  The  cross-wind  model  is  presented  in  Appendix  C 
and  is  summarized  below.  Justification,  if  any,  for  the  model 
awaits  the  collection  of  wind  data  (Volume  TT  of  this  report). 
Existing  data  (discussed  in  Appendix  C)  do  support  the  model  even 
though  the  derivation  of  the  model  contains  some  inconsistencies. 

The  cross-wind  model  was  formulated  using  winds  measured  by 
aerovanes  mounted  on  a 1500- foot  tower  (Ref.  31).  The  cross- 
wind  p.d.f.,  Pc(wl.),  in  terms  of  the  mean  wind  speed  aloft  u may 
be  written  (Appendix  C) : 

■)  7 

n r,,  i - 1 *irw  “/4u 
Pc'.wc-'  u e c 
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for  wc  > 0.  Using  p=18.6  knots  (one  of  the  observed  values  in 
Ref.  31), 


t \ n nr?o  “0.00227w  { \ 

pc(wc)  = 0.0538e  c . (29) 

The  VAS  measures  wind  velocity  50  feet  above  the  ground  at 
the  approach  end  of  the  active  runway.  The  algorithm  for  deter- 
mining "vortex  safe"  or  "green  light"  conditions  in  the  MM  to  run- 
way threshold  region  relies  on  the  one-minute  average  wind  vector 
observation  lying  outside  an  ellipse  (Tig.  5)  that  is  aligned  with 
the  runway  coordinate  system.  The  calculation  of  surface  cross- 
wind  probability  density  is  straightforward,  involving  the  develop- 
ment of  a conditional  surface  wind  rose  (based  on  the  wind  vector 
not  being  inside  the  VAS  ellipse),  and  then,  performing  the  vari- 
able transformation  on  the  conditional  surface  wind  rose. 

The  procedure  is  less  clear  for  inferring  a cross-wind  p.d.f. 
aloft  given  a "red  light/green  light"  observation  of  surface  winds. 
Little  information  exists  that  can  be  brought  to  bear  on  this 
problem.  Generalizations  can  be  made:  steering  or  veering  of  the 
wind  occurs  with  increasing  altitude,  and  wind  speed  increases  with 
altitude.  Specific  data  that  could  be  used  to  define  an  analog  to 
the  surface  ellipse  for  use  at  altitude  do  not  appear  to  exist. 

In  the  absence  of  such  specific  detail,  a simple  statement  is 
postulated:  if  the  VAS  registers  a "green  light,"  the  wrinds  aloft 

are  also  outside  of  the  VAS  ellipse.  There  are  several  questions 
and  additional  assumptions  implicit  in  the  use  of  this  postulate 
(discussed  at  length  in  Appendix  C) . However,  limited  data  col- 
lected at  O'Hare  appear  to  be  in  agreement  with  the  postulate. 

The  verification  of  this  correspondence  is  one  of  the  subjects  of 
Volume  II.  Figure  13  shows  the  two  p.d.f. 's:  p (w  ) (equation 
(29))  for  the  cross-wind  magnitude  with  no  wind  information  from 
the  VAS,  and  p^fw^/GLa)  for  the  cross-wind  magnitude  when  the  VAS 
registers  a "green  light." 


PROBABILITY  DENSITY 


4.5  analytic:  modi:! 


Sections  4.1  through  4.4  describe  an  analytic  model  for  cal- 
culating the  probability  of  a hazardous- vortex  encounter.  A com- 
puter routine  for  doing  the  calculations  is  given  in  Appendix  D. 
The  steps  in  the  calculations  are  summarized  below. 

The  initial  vortex  strength  is  calculated  from: 


ro  = mbe+r int ’ 


where  b£  is  the  wingspan  of  the  vortex-encountering  aircraft, 
decay  of  the  vortex  strength  is  given  by  (equation  (8)) 


The 


r = 9-58 
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for  d > 47.9  b . The  hazard  radius  is  then  calculated: 


R = 
o 


2irfpU 


If  Ro>  bc/2,  then  a hazardous-vortex  encounter  is  possible  depend- 


ing on  the  location  of  the  following  aircraft  with  respect  to  the 
vortices.  The  horizontal  encounter  probability  is  (equation  (20)) 
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the  vertical  encounter  probability  is  (equation  (21)) 
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so  that  the  probability  that  the  vortex- encounter ing  aircraft  is 
within  the  hazardous- vortex  encounter  cross  section  is: 


PHE  PVE  dw’ 


where  p (w)  is  the  p.d.f.  for  the  wind.  The  probability  of  a 
hazardous- vortex  encounter  is: 


P,  = P,  P 
he  he 


where  (equation  (14)): 
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5,  RESULTS  OF  SAFETY  ANALYSIS  FOR  ALL  AIRCRAFT 


In  this  section,  the  results  of  the  V'AS  safety  analysis  are 
presented.  The  analytical  process  consists  of  seven  fundamental 
steps  : 

1.  Calculation  of  the  ordinary  risks  accepted  in  the  present 
system  using  3/4/5-nautical-mile  spacings. 

2.  Establishment  of  a baseline  probability  as  a maximum 
acceptable  risk. 

3.  Calculation  of  the  new  risks  to  be  expected  when  using 
3-nautical-mile  spacings  during  VAS  green-light  conditions. 

4.  Formulation  of  operational  guidelines  to  maintain  the 
risks  with  VAS-reduced  spacings  at  or  below  the  baseline  risk. 

5.  Sensitivity  analysis  on  the  conclusions  to  determine  if 
the  conclusions  are  adversely  affected  by  potential  errors  in  the 
models  and/or  parameters. 

6.  Comparison  of  the  conclusions  to  the  available  body  of 
physical  evidence  of  risk  and  acceptability. 

7.  Formulation  of  a refined  set  of  operational  guidelines 
and  recommended  actions. 

The  model  presented  in  Section  4 is  used  for  the  calculations.  The 
fraction  f of  the  roll-control  capability  (Section  4.1.2)  is  taken 
as  0.378. 

5.1  RISKS  IN  PRESENT  SYSTEM 

In  any  probability  analysis  dealing  with  small  numbers,  it  is 
important  to  interpret  the  results  in  a meaningful  way.  Because  of 
the  many  conservative  assumptions  enacted  in  defining  the  model  and 
its  parameters,  the  absolute  probabilities  themselves  have  limited 
meaning.  Comparison  among  the  various  probab i 1 i t ies  are  expected  to 
be  meaningful  however  (possible  exception  is  discussed  in  Sec.  5.2). 

The  largest  probability  for  a hazardous- vortex  encounter  at 
the  current  separations  occurs  for  the  PA- 28  following  3 nautical 
miles  behind  the  PC  - 8 . Table  4 shows  that  the  hazard  radius  shrinks 
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to  zero  at  a separation  of  6 nautical  miles.  (Note  that  this  dis- 
cussion refers  to  the  region  near  the  OM  where  the  separation 
standard  is  3 nautical  miles;  the  separation  standard  calls  for 
4 nautical  miles  near  the  MM  due  to  the  possibility  of  a vortex 
stalling  in  ground  effect.)  All  probabilities  will  be  compared 
with  the  DC-8/PA-28  probability  of  a hazardous -vortex  encounter. 

Table  6 shows  the  risks  at  the  OM  for  the  various  aircraft 
pairs  at  the  current  3/4/5  separations  relative  to  the  risk  for 
the  DC-8/PA-28  pair.  Most  of  the  entries  are  zero.  According  to 
the  model,  the  risk  of  a hazardous  - vortex  encounter  is  zero  for  the 
airliners  following  other  airliners  at  the  current  separations  and 
nonzero  only  for  the  Small  category  of  aircraft  following  certain 
airliners.  However,  no  accident  is  known  to  have  been  caused  by 
vortices  when  any  of  the  respective  aircraft  were  separated  by  the 
appropriate  3,  4 or  5 nautical  miles. 

5.2  BASELINE  PROBABILITY 

The  situations  represented  by  the  nonzero  probabilities  in 
Table  6 (the  PA-28  5 nautical  miles  behind  a B-747,  DC-10,  L-1011, 
DC-8H,  or  B-707II;  the  Lcarjet  5 nautical  miles  behind  a B-747;  the 
PA-  28  3 nautical  miles  behind  a DC-8,  B-  707,  or  B - 7 2 7 ; and  the 
Learjet  3 nautical  miles  behind  a DC-8  or  B-707)  occur  frequently 
in  the  vicinity  of  the  OM.  Since  these  situations  occur  frequently 
(particularly  in  VFR  conditions)  and  with  nonzero  probability  of  a 
vortex  hazard,  these  risks  are  considered  to  be  acceptable  as  no 
hazardous  - vortex  encounters  are  known  to  have  occurred.  Refer- 
ence 1 (and  Section  2)  examined  these  situations  in  great  detail, 
and  found  no  apparent  hazard. 

The  baseline  probability,  N,  is  defined  as  the  largest  of  the 
nonzero  probabilities  (DC-8/PA- 28) . If  the  probability  of  a vortex 
encounter  for  any  leader/follower  pair  is  found  to  be  less  than  or 
equal  to  N,  the  situation  will  be  assumed  to  be  safe.  If  the 
probability  of  a hazardous -vortex  encounter  were  greater  than  N, 
it  would  not  necessarily  mean  that  a hazardous  condition  existed- - 
it  would  mean,  however,  that  a vortex  encounter  might  occur,  and 


the  degree  of  hazard  could  depend  on  how  much  larger  the  probabil- 
ity was  compared  with  N (and  to  the  actual  roll  recovery  capability 
of  the  following  aircraft).  Also,  it  does  not  mean  that  1 in  li  ' 
landings  of  a I’A-28  behind  a DC-8  will  result  in  a hazardous-vor- 
tex  encounter;  the  model  and  its  parameters  have  been  purposely 
chosen  to  be  conservative. 

There  is,  however,  a possible  problem  with  comparing  two  con- 
servative estimates.  Absolute  conservatism  of  estimates  of  two 
elements  does  not  necessarily  imply  conservation  in  their  compari- 
sons. For  example,  see  Fig.  14.  Consider  a proposed  case  where 
reduction  to  3 nautical  niles  under  VAS-green  conditions  is  recom- 
mended. Both  the  estimated  baseline  and  the  estimated  proposed 
probabilities  of  a hazardous  vortex  encounter  are  conservative 
estimates  of  the  actual  baseline  and  proposed  probabilities.  Based 
on  these  estimates  one  would  claim  that  the  hazard  associated  with 
the  proposed  case  is  below  that  associated  with  the  baseline  case 
(because  proposed/estimated  < baseline/estimated)  and  hence  the 
proposed  case  is  acceptable.  However,  it  is  possible  that  the 
actual  probability  of  the  proposed  case  is  greater  than  the  actual 
probability  of  the  baseline  case  (as  is  shown  in  Fig.  14),  and 
hence  the  situation  represented  by  the  proposed  case  is  more 
hazardous  than  the  baseline  case.  To  determine  the  consequences 
(if  any)  of  the  situation  just  described  requires  detailed  analysis 
and  is  deferred  to  Volume  III. 

If  the  reader  is  concerned  about  the  possible  problem  and  its 
consequences,  the  analysis  in  Section  6 avoids  the  problem  by  in- 
troducing a VAS  which  excludes  reduced  separations  for  the  Small- 
category  aircraft.  The  remainder  of  Section  5 assumes  that  the 
problem  is  resolved  in  favor  of  the  analytical  approach  used  here- 
in. Even  if  the  problem  of  comparing  conservative  estimates  does 
apply  in  this  analysis,  most  of  the  results  in  this  section  will 
not  be  altered. 

5.3  RISKS  WITH  VAS 

The  probabilities  of  a hazardous-vortex  encounter  are  shown 
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in  Table  7 for  VAS-green  conditions.  Here,  all  the  aircraft  are 
spaced  at  3 nautical  miles,  and  the  probabilities  are  given  rela- 
tive to  N,  the  baseline  value.  Only  one  entry  exceeds  the  baseline 
probability;  a value  of  1 . 8N  is  found  for  the  PA-28,  3 nautical 
miles  behind  the  B-747.  This  probability  reduces  to  1 . ON  for  a 
separation  of  3.5  nautical  miles.  Figure  15  shows  how  the  proba- 
bility (in  units  of  N)  varies  as  a function  of  the  distance  of  the 
PA- 28  from  the  runway  threshold,  given  that  the  PA  - 2 8 is  3 nautical 
miles  behind  the  B-747  in  VAS-green  conditions.  If  the  OM  is  more 
than  29,000  feet  (4.7  nautical  miles)  from  the  runway  threshold, 
then  the  probabilities  can  exceed  the  baseline  value. 

It  is  important  to  reiterate  that  decreased  separations  will 
be  permitted  only  when  the  VAS  displays  a green  light.  During 
VAS-red  conditions,  the  model  indicates  probabilities  of  a hazard- 
ous vortex  encounter  in  excess  of  N.  Because  of  the  difference  in 
flight  speeds,  it  is  unlikely  that  many  of  the  aircraft  pairs  can 
maintain  a 3-nautical-mile  separation  (e.g.,  B-  74 7/PA- 28 ) . The 
uniform  3- naut ical -mi le  separation  is  intended  only  inside  the  OM. 
In  most  situations  the  aircraft  are  in  trail  with  the  present 
standards  well  before  reaching  the  OM  so  the  minimum  separations 
obtainable  using  VAS  will  be  approximately  3,  3^,  and  \h  nautical 
miles. 

5.4  OPT RATIONAL  GUIDELINES 

To  maintain  the  risks  with  VAS-reduced  spacings  at  or  below 
the  baseline  risk  N,  certain  operational  guidelines  are  required. 
Table  7 shows  that  nonzero  probabilities  are  obtained  for  the  Small 
category  of  aircraft  following  most  of  the  airliners.  Reduced 
separations  are  permissible  during  VAS-green  conditions  as  the 
vortices  of  the  preceding  aircraft  are  descending  and/or  translat- 
ing out  of  the  path  of  the  following  aircraft.  There  may  be  prob- 
lems if  these  Small  aircraft  (as  well  as  the  R-727,  DC-9,  and 
B-737  behind  a B-747)  are  turned  on  at,  or  allowed  to  descend  to, 
lower  altitudes  than  the  leading  aircraft.  During  the  trails  of 
the  VAS,  the  Small  category  of  aircraft  should  be  positioned  at  3 
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FIGURE  15.  NORMALIZED  PROBABILITY  OF  VORTEX  ENCOUNTER  FOR  A PA-28 
THREE  NAUTICAL  MILES  BEHIND  A B-747  UNDER  VAS-GREEN  CONDITIONS  AS 
A FUNCTION  OF  THE  DISTANCE  OF  THE  PA- 28  FROM  THE  RUNWAY  THRESHOLD 


nautical  miles  behind  other  categories  only  if  both  aircraft  are 
turned  on  the  ILS  at  or  beyond  the  OM,  and  both  are  making  preci- 
sion (i.e.,  localizer  and  glide  slope)  approaches.  Because  of  the 
generally  slower  speeds  of  the  Small  aircraft,  prohibiting  short 
finals  will  cause  the  3-naut ical -mi  1 e separation  at  the  OM  to  in- 
crease gradually  to  4 nautical  miles  at  the  MM. 

Although  the  probability  of  a vortex  encounter  for  the  B-747/ 
PA-28  pair  at  3 nautical  miles  is  found  to  be  in  excess  of  the 
baseline  probability  N,  the  conservative  assumptions  throughout 
the  analysis  should  preclude  the  necessity  for  segregating  the 
PA-28  (and  other  civil  single-engine  or  light-twin  aircraft)  from 
following  the  B-747  at  3 nautical  miles.  The  problem  (if,  indeed, 
there  is  a problem)  exists  only  in  the  vicinity  of  the  OM  where 
the  PA- 28  will  have  ample  room  to  maneuver. 

5.5  SENSITIVITY  ANALYSIS 

A sensitivity  analysis  has  been  done  to  determine  the  effect 
on  the  normalized  encounter  probability  and  f of  possible  errors 
in  the  models  and/or  parameters.  Details  are  presented  in  Appendix 
E.  Although  the  probability  values  do  change  depending  on  the 
details  of  the  wind,  navigation,  vortex-descent,  vortex  decay,  and 
rol I- recovery  capability  models,  the  results  and  conclusions  con- 
cerning the  safety  of  the  VAS , however,  do  not  appear  to  require 
alteration.  All  the  calculations  were  done  with  the  following  air- 
craft at  the  OM  as  this  is  the  location  at  which  the  encounter 
probability  is  maximized.  The  OM  was  assumed  to  be  7 nautical 
miles  from  the  runway  as  this  also  maximizes  the  probability. 

5.6  OTHER  EVIDENCE 

Other  evidence  was  examined  to  look  for  support  or  contradic- 
tion of  the  conclusions  of  the  safety  analysis.  Three  independent 
observations  lend  credence  to  the  analytic  results:  wake  vortex 
incident  reports  analyzed  by  the  British  Civil  Aviation  Authority, 
accident  data,  and  VFR  spacings. 

The  Civil  Aviation  Authority  of  Great  Britain  has  been 


analyzing  wake  vortex  incidents  (no  accidents)  reported  in  the 
United  kingdom  since  1972  (Refs,  b-8  and  Section  3.1).  Through 
September  1977,  353  incidents  have  been  tabulated.  When  the  docu- 
mented (73)  incidents  occurring  between  runway  threshold  and  an 
altitude  of  1000  feet  above  ground  level  were  correlated  with  the 
wind  measured  at  a height  of  30  feet  (the  centerfield  sensor  at 
the  London  Heathrow  International  Airport),  only  3 incidents  were 
found  when  the  winds  exceeded  the  VAS  wind  criterion  in  Figure  5 
(Refs.  7 and  8).  All  severe  incidents  (reported  bank  angles  in 
excess  of  30  degrees)  occurred  in  wind  conditions  well  inside  the 
ellipse  (VAS-red  conditions).  The  three  incidents  which  occurred 
in  wind  conditions  outside  the  ellipse  (VAS-green  conditions)  were 
within  2 knots  of  the  wind  criterion.  The  2 incidents  farthest 
outside  the  ellipse  occurred  at  an  altitude  of  1000  feet  (one  had 
a reported  bank  angle  between  10  and  30  degrees,  and  the  other  less 
than  10  degrees).  The  third  incident  occurred  with  centerfield 
winds  just  outside  the  ellipse  (less  than  0.25  knots  from  the 
boundary)  and  reported  a roll  excursion  of  less  than  10  degrees. 

The  Heathrow  centerfield  winds  were  averaged  over  a 10-minute 
period  and  the  results  do  not  contradict  the  VAS  concept  and  the 
safety  analysis.  The  VAS,  however,  uses  1 -mi note-averaged  winds. 

An  analysis  of  vortex-caused  accidents  (Ref.  1 and  Section 
2.2.1)  noted  that  those  accidents  which  occurred  between  the  MM 
and  DM  have  distinct  characteristics.  The  vortex  encounters  seemed 
to  be  random,  and  occurred  only  because  the  Small  aircraft  were 
approaching  at  a low  altitude  with  poor  lateral  navigation.  These 
aircraft  apparently  would  not  have  encountered  vortices  if  the 
aircraft  were  on  a precision  approach. 

The  chance  of  a vortex  encounter  increases  as  interarrival 
aircraft  separations  decrease.  VFR  operations  oftentimes  are  con- 
ducted with  interarrival  separations  which  are  less  than  the  1 FR 
minimum  separations.  But  these  smaller  spacings  are  flown  safely 
in  VFR.  Of  course,  VFR  and  IFR  operations  are  not  the  same;  in 
VFR  the  pilot  accepts  more  risk  as  he  is  more  aware  of  his 
environment.  The  VAS  proposes  a 3-nautical-mile  minimum  and  only 
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during  specific  wind  conditions. 


5.7  CONCLUSIONS  AND  RECOMMENDED  ACTIONS 

The  VAS  as  described  herein  is  intended  for  independent 
s ingle  - runway  landings.  Closely  spaced  parallel  runways  and  runway 
intersections  will  require  a case-by-case  examination  to  determine 
how  VAS  should  operate.  For  example,  an  8-knot  cross  wind  indi- 
cates VAS-green  conditions,  and  separations  may  be  reduced  to  a 
uniform  3 nautical  miles.  However,  if  the  vortices  are  blown 
toward  a closely  spaced  parallel  runway,  the  wind  criterion  that 
insures  vortices  cannot  pose  a threat  to  an  aircraft  landing  on 
the  same  runway  may  not  insure  that  the  threat  is  nonexistent  for 
an  aircraft  landing  on  a downwind  parallel  runway.  Either  the 
wind  criterion  will  need  to  be  modified  or  operational  restrictions 
will  need  to  be  enforced  to  prevent  a hazardous  condition. 

The  analysis  has  not  considered  other  than  straight-in  ap- 
proaches that  might  occur  under  IFR  (e.g.,  sidestep  or  circling 
approaches).  These  other  approaches  are  not  necessarily  unsafe 
with  the  VAS;  indeed,  such  extensions  need  to  be  done  at  some 
future  time,  but  the  analysis  herein  addressed  only  straight-in 
IFR  approaches. 

The  results  of  the  analysis  show  that,  subject  to  verifica- 
tion of  several  assumptions,  uniform  3-nautical-mile  separations 
between  the  runway  threshold  and  the  outer  marker  may  be  used 
safely  when:  (1)  the  VAS  registers  a green  light,  and  (2)  all 
aircraft  are  restricted  to  an  ILS  approach  between  the  OM  and 
touch  down.  Short  finals  and  VOR  or  localizer  approaches  with 
reduced  separations  should  not  be  permitted  (short  finals  are  pre- 
cluded in  IFR);  the  flight  paths  of  these  non-precision  approaches 
act  to  increase  the  probability  of  a vortex  hazard.  The  normal 
vortex  descent  and  translation  with  the  wind  which  decrease  the 
hazard  probability  for  precision  approaches  can  move  the  vortices 
into  the  path  of  an  aircraft  on  a short  final  or  on  a VOR  or 
localizer  approach.  Extensive  data  have  demonstrated  the  safety 
in  the  threshold  to  middle  marker  region;  probability  analysis 
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based  on  state-of-the-art  knowledge  of  vortex  behavior  extends  the 
safe  region  from  middle  to  outer  marker.  For  many  of  the  aircraft 
pairs  a simplified  analysis  is  appropriate  that  avoids  the  use  of 
unverified  assumptions  (see  Section  6). 

The  results  herein  are  derived  from  an  analytical  model.  A 
data-col lection  program  is  underway  to  verify  the  mathematical 
analysis.  Appendix  F describes  the  data-col lect ion  effort  at 
O'Hare.  Until  the  confirmation  tests  are  completed,  the  analysis 
is  only  as  good  as  the  model.  However,  the  model  is  conservative 
and  has  been  shown  to  be  relatively  insensitive  to  the  assumed 
parameters . 

When  the  YAS  indicates  reduced  separations  may  be  used,  the 
state  of  knowledge  of  vortex  behavior  is  consistent  with  safely 
applying  the  reduced  separations  out  to  the  outer  marker.  If  the 
tests  at  O'Hare  confirm  the  model  predictions,  the  FAA  should  be 
justified  in  commissioning  the  VAS  as  a tool  for  decreasing  delays 
at  the  major  airports. 


6,  SAFETY  ANALYSIS  FOR  HEAVY/HEAVY  AND  HEAVY/LARGE 


The  two  main  unresolved  issues  of  the  approach  presented  in 
Sections  4 and  5 were  the  wind  model  (Section  4.4  and  Appendix  C) 
and  the  potential  problem  related  to  the  comparison  of  two  conserva- 
tive estimates  (Section  5.2).  If,  however,  the  VAS  were  to  be  used 
to  reduce  only  the  Heavy/Heavy  and  Heavv/Large  standards  (omitting 
the  Heavy/Small  and  Large/Small  cases),  then  the  two  unresolved 
issues  will  not  be  considered  in  the  safety  analysis  as  will  be 
shown  in  this  section.  The  discussion  below  addresses  the  safety 
of  3-nautical-mile  separations  for  jetliners  during  VAS-green 
conditions,  and  then  expands  the  analysis  to  include  all  Large 
aircraft . 

6.1  HAZARD  MODEL 

In  Section  4.1.1,  it  was  argued  that  a potentially  hazardous- 
vortex  encounter  can  occur  only  if  _>  b^/2  and  if  a following 
aircraft  were  fully  located  within  the  circle  defined  by  the  hazard 
radius  Rn-  The  probability  of  a hazard  is  defined  as  zero  when 


If  one  accepts  this  definition,  then  for  all  those  pairs  of  air- 


craft with  R < 
o 


b /2  at  3 nautical  miles,  the  standards  can  be 
e 


reduced.  The  two  unresolved  issues  do  not  apply  to  these  cases. 

Using  f = 0.378  (one-half  the  maximum  roll  authority  required 

safely  to  fly  the  current  separation  standards  at  the  assumed 

speeds)  and  equation  (14)  (or  equation  (12)  if  necessary)  , the 

distances  at  which  R becomes  smaller  than  b /2  are  calculated  and 

o e 

presented  in  Table  4.  Most  of  the  calculated  distances  are  less 
than  3 nautical  miles;  for  these  aircraft  pairs,  the  model  in- 
dicates that  the  separation  standards  can  be  set  at  3 nautical 
miles  between  the  outer  and  middle  marker  with  or  without  the  VAS. 
The  VAS  will  be  required  to  indicate  when  3 nautical  miles  may  be 
used  between  the  middle  marker  and  the  runway  threshold  however. 
The  Smal 1 -category  aircraft  following  Heavy-anJ  Large- category 
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aircraft  and  the  B-727,  B-737,  and  DC-9  following  the  B-747  have 

yielded  calculated  distances  in  excess  of  3 nautical  miles.  It 

will  be  shown  that  the  condition  of  R < b /2  can  be  expanded  to 

o e 

embrace  the  B-727,  B-737,  and  DC-9  following  the  B-747,  and  that 
the  VAS  can  be  applied  to  the  operations  of  Heavy/Heavy  and  Heavy/ 
Large  combinations  of  aircraft. 

Recalling  equations  (5)  and  (8),  the  hazard  radius  is 

R = r 
0 2nfp  II 
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for  dC. /b  A > 9.58. 
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Introducing  equation  (13)  and  rearranging, 
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The  following  relationships  can  be  inferred  from  the  above 
equat ion : 

a.  Rq  increases  if  f,  p,  or  U decreases. 

b.  Rq  increases  if  A/C.  or  the  estimate  of  vortex  decay 
(indicated  by  9.58,  the  breakpoint  of  the  decay  curve-- 
Fig.  6)  increases. 

For  a given  aircraft  pair,  b , b , and  d are  constants.  Since  m 

g ^ 

and  T-  are  selected  as  ultimate  or  maximum  values  (Appendix  A), 
in  i a 

they  may  be  considered  constants.  A p of  0.06  (airliners)  or  0.08 
(general-aviation  aircraft)  is  a lower  bound  (Section  4.1.1),  so 
that  it  also  may  be  considered  a constant.  The  9.58  value  is  an 
upper  limit  on  the  decay  parameter  (Fig.  6)  except,  possibly,  for 
the  B-727.  For  the  calculations  to  follow,  12.0  is  used  for  the 
B-727  decay  parameter  (the  suspect  B-727  data  in  Fig.  6 are  in- 
cluded), and  9.58  is  used  for  all  other  vortex -generat ing  or  lead 
aircraft . 
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Thus  the  vortex  decay  parameter  may  also  be  considered  a constant. 


Setting  equal  to  bp/2, 

. tA/Cd  K * rintK  (9-581  . 

11  irp  db^ 

Table  3 is  obtained  by  setting  d equal  to  the  current  separation 
standards;  Table  8 is  obtained  by  setting  d equal  to  3 nautical 
miles--many  of  the  entries  are  the  same  as  in  Table  3. 

The  largest  f value  in  Table  8 is  0.458  for  the  B-747/DC-9 
pair.  If  this  f value  is  acceptable  (i.e.,  assuming  the  aircraft 
ever  finds  itself  in  a vortex,  it  will  require  no  more  than  46 
percent  of  an  aircraft's  roll-control  authority  to  counteract  the 
rolling  motion  of  a vortex)  , then  the  3-nautical-mile  standard 
can  be  applied  under  the  assumption  that  Rq  < b /2  implies  no 
vortex  hazard.  It  is  noted  (see  Table  3)  that,  according  to  the 
hazard  model  f values  greater  than  0.458  are  required  now  at  the 
current  separation  standards  and  assumed  speeds  for  the  DC-8/Lear- 
jet  and  for  the  general-aviation-type  aircraft  following  Heavy  and 
Large  aircraft. 

6.2  SENSITIVITY  ANALYSIS 

The  purpose  of  this  section  is  to  examine  the  sensitivity  of 
the  parameters  in  the  expression 

= CA/Cl)  (mb,  + rintJhg  (9.58)  . 

U 71  p 

As  indicated  for  R in  Section  6.1,  the  elements  of  importance  are 
(A/Cl)  and  U.  Differentiating 

4f  ■ 3TA7qr  k(A/c0  * ro  4U 


L 
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d.2.1  B - 74  7/1)0  9 

Consider  first  the  B-747/DC-9  combination  as  it  gives  the 
largest  f value  for  the  airliners  at  3-naut ical -mi le  separation. 

The  maximum  non-dimensional  roll  rate,  p,  is  assumed  to  be  0.06 
as  this  is  the  minimum  value  for  airliners.  However,  the  DC-9  has 
considerably  more  lateral  roll  power.  Reference  32  indicates  that 
in  landing  configuration  with  rudder-fixed  p varies  from  about  0.07 
at  100  knots  to  0.25  at  180  knots;  with  the  outboard  spoiler  in- 
operative the  value  at  180  knots  falls  to  0.16.  At  the  speed 
assumed  in  the  calculations  for  Table  8,  p is  about  0.12.  Thus,  a 
more  realistic  value  of  f will  be  (0 . 4 58 ) (0 . 06) / (0 . 12 ) = 0.229, 
and  the  B-747/DC-9  combination  no  longer  represents  the  worst-case 
situation  (even  with  different  speeds). 


6.2.2  B-747/B-757 

The  B-747/B-737  pairing  has  yielded  the  second  largest  f value 
in  Table  8.  Aerodynamic  coefficients  in  Ref.  33  indicate  that  p 
is  about  0.10  for  the  B-737.  Thus,  a more  realistic  value  of  f is 
0.265,  and  now  the  B-747/B-737  combination  no  longer  represents 
a worst-case  situation. 

6.2.3  B- 747/B-727 

Since  the  B-747/DC-9  and  B-747/B-737  pairings  have  been 
eliminated  as  contenders  for  the  worst-case  situation,  the  B-747/ 
B-727  is  next  in  line.  The  pertinent  aerodynamic  coefficients  for 
determining  p for  the  B-727  could  not  be  obtained.  Thus,  f = 0.400 
is  taken  to  be  the  largest  f value  required  by  an  airliner  to  fly 
safely  3-naut ica 1 -mi lc  spacings  between  the  outer  and  middle 
marker.  This  figure  is  less  than  some  f values  calculated  using 
today's  separation  standards  (see  Table  3). 

An  upper  bound  on  the  uncertainty  in  f is  found  by  evaluating 
equation  (30)  using  the  largest  positive  deviation  in  A/C.  and  the 
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largest  negative  deviation  in  l).  A/Cj  is  approximated  by  5.0  in 
all  calculations.  A more  precise  value  for  the  B-747  is  5.1  (Cj 
is  about  1.4  for  a landing  B-747  (Ref.  21)  and  A is  the  wingspan 
divided  by  the  mean  aerodynamic  chord,  195.7/27.3  = 7.2).  The 
approach  speed  for  the  B-727  varies  between  about  105  and  130  knots 
for  landing  weights  from  operational  empty  weight  plus  reserves  to 
maximum  landing  weight.  Using  the  slowest  speed, 


Af  = jp™  f5-1  ' 5-0) 

= 0.008  + 0.056 
= 0.064. 


0.400 

205.8 


(177.2 


205.8) 


The  maximum  f required  for  3-nautical-mile  separations  is  f + Af 
or  0.464,  less  than  half  the  available  roll-control  power.  However, 
f values  greater  than  0.464  are  found  for  some  aircraft  pairs  using 
today's  standards  (Table  3).  This  suggests  that  VAS-reduced 
separations  for  jetliners  are  safe. 

6.2.4  Other  Aircraft  Pairs 

Similarly,  the  remaining  leader/follower  combinations  can  be 
examined.  However,  the  B-747/B-727  case  discussed  above  leads  to 
the  largest  value  of  f + Af. 

The  argument  presented  herein  cannot  he  used  with  following 
Small  aircraft.  The  PA-28  would  require  f values  in  excess  of  1.0 
following  the  B-747,  DC-8,  and  B-707  at  3 nautical  miles  as  may  be 
seen  by  appropriately  scaling  the  entries  in  Table  3. 

6.3  VAS  FOR  I1LAVY  AND  LARGE  AIRCRAFT 

In  Section  6.2,  it  is  shown  that  the  decay  of  vortices  alone 
will  permit  3-nautical-mile  separations  between  the  outer  and  mid- 
dle markers  for  the  airliners  or  transport-category  jet  aircraft. 
Using  VAS-reduced  spacings  for  a Heavy  or  a Large  airliner  follow- 
ing a Heavy  probably  can  be  implemented,  but  effectively  introduces 
a new  category  for  air  traffic  controllers  to  contend  with  (Heavy, 
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Large,  ana  Small  versus  Heavy,  Large  jetliners,  other  Large,  and 
Small).  Although  the  "other  Large"  category  aircraft  operate  less 
frequently  into  the  large  airports,  it  will  obviously  be  easier 
for  the  controllers  if  all  Large  aircraft  use  VAS-reduced  spacings. 
In  this  section,  it  is  shown  that  VAS-reduced  separations  may  be 
used  safely  by  lleavy/Heavy  and  Heavy/Large  combinations;  the 
analysis  does  not  depend  on  the  unresolved  issues  in  Section  5; 
namely,  the  wind  model  and  the  possible  problem  with  comparing 
conservative  probability  estimates. 

Given  an  acceptable  f value,  the  minimum  approach  speeds  con- 
sistent with  a safe  vortex  encounter  (P^  = 0)  can  be  calculated 
for  the  various  Large  aircraft.  However,  what  is  an  acceptable 
value  of  f?  Table  3 shows  that  the  maximum  f required  for  the 
current  separation  standards  is  0.756  for  the  PA-28  3 nautical 
miles  behind  a DC-8.  To  be  conservative,  a slow  speed  is  used  in 
those  calculations  for  the  PA-28;  a more  typical  PA-28  speed  at 
the  outer  marker  is  about  95  knots.  Using  95  knots,  f reduces  to 
0.519.  For  the  purposes  of  this  analysis,  0.519  will  be  used  for 
today's  separations,  rather  than  the  more  conservative  0.758.  Using 

f . (A/CU  (mbe  * Unt)  \ (9'58>  , 

^ up  dh^ 

a parametric  equation  can  be  found  for  determining  the  minimum 
flight  speed  U required  to  give  f = 0.519  for  the  various  Large 
aircraft.  The  worst-case  situation  is  when  a B-747  is  the  vortex- 
generating or  lead  aircraft.  Selecting  the  B-747  specifies  m, 

^int’  '3g’  an^  A/Cl*  With  some  rearranging  of  terms,  the  above 
equation  becomes 

(5.1) (9.58) (195.7)  119.56  b + 1148.6) 

U = Mp(3) (60.76) (0.519)be 

( 1 9 . 5 6 bc  + 11 4 8. 6 ) 


0.322 


Figure  lb  shows  U versus  b for  f = 0.519.  The  two  curves  are 

/s  C 

for  p = 0.06  and  0.08.  For  a specific  aircarft  (e.g.,  b^  and  p) , 
the  curves  show  the  minimum  speed  required  to  give  f = 0.519;  higher 
speeds  lead  to  lower  values  of  f.  For  those  aircraft  for  which 
aerodynamic  coefficients  arc  available,  the  X's  mark  the  calculated 
minimum  U.  For  these  aircraft,  as  long  as  the  speed  U is  greater 
than  or  equal  to  the  plotted  value,  f will  be  less  than  or  equal 
to  0.519.  The  interpretation  of  Fig.  16  is  as  follows.  The  model 
(Section  4.1.1  and  4.1.2)  gives  a zero  probability  of  a hazardous- 
vortex  encounter  at  VAS-reduced  spacings  and  an  assumed  safe  f = 
0.519  if,  for  instance,  the  Learjet-25  passes  the  outer  marker  with 
a speed  of  at  least  90  knots  and  the  DC-9  at  least  50  knots  (since 
DC-9s  require  a speed  in  excess  of  50  knots  to  fly,  they  certainly 
pass  the  outer  marker  with  a speed  of  at  least  50  knots). 

For  those  aircraft  where  p is  not  available,  the  minimum 
flight  speed  at  the  outer  marker  required  to  give  = 0 may  be 
read  from  the  appropriate  curves  in  Fig.  16.  excluding  military 
aircraft  (fighters  are  designed  to  achieve  at  least  p = 0.09),  all 
the  Large  aircraft  can  be  checked  to  find  the  largest  minimum  U. 

For  example,  after  the  Learjet  models,  the  next  larger  wingspans 
are  42.3  feet  for  the  Aerospatiale  Corvette  and  42.9  feet  for  the 
Dassault  Falcon-10.  These  aircraft  have  a p of  at  least  0.08,  anu 
thus,  lead  to  a minimum  U of  111  and  110  knots,  respectively. 
Similarly,  checking  each  of  the  Large  aircraft  and  using  the  ap- 
propriate p leads  to  a largest  minimum  U of  111  knots. 

The  calculations  above  have  shown  that  VAS-reduced  spacings 
are  safe  between  the  outer  and  middle  marker  for  the  airliners  and 
for  the  non-airliner  Large  aircraft  if  their  speeds  at  the  OM  are 
at  least  111  knots.  Since  approach  speeds  are  usually  in  excess  of 
120  knots  at  the  outer  marker,  the  use  of  VAS-reduced  spacings  may 
be  applied  for  Heavy/Heavy  and  Heavy/Large  combinations. 

If  a maximum  I of  0.519  is  acceptable  (recall  that  some  of 
today’s  standards  reflect  larger  f values-  Table  3),  the  VAS  can 
be  used  to  permit  3- naut ical-mi le  separations  in  VAS-green  condi- 
tions for  Heavy  and  Large  aircraft  following  a Heavy.  It  must  be 
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noted  that  a maximum  f of  0.519  does  not  mean  for  instance,  that 
every  aircraft  following  3 nautical  miles  behind  a B-747  will  re- 
quire the  use  of  half  its  roll-control  power  to  conduct  its 
approach.  According  to  the  model,  if  an  aircraft  ever  finds  itself 

in  the  vortex  of  a B-747  3 nautical  miles  ahead,  the  extreme  upper 
theoretical  limit  on  the  f required  will  be  about  0.5,  within  the 
control  capability  of  aircraft.  Pending  further  analyses  and  data 
collection  (Volumes  III  and  II,  respectively),  the  Smal 1 -category 
aircraft,  if  this  restricted  VAS  is  adopted,  retain  their  present 
separation  standards  regardless  of  the  VAS  condition. 
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7,  SUMMARY 


The  discussions  in  Sections  5 and  6 have  presented  al- 
ternative safety  analyses  of  the  VAS.  Both  analyses  used  the  same 
model  of  a vortex  hazard  (Sections  4.1.1  and  4.1.2).  At  this 
point,  the  two  analyses  follow  different  arguments.  The  restricted 
VAS  in  Section  6 notes  that  vortex  decay  is  the  primary  mechanism 
outside  the  middle  marker  for  avoiding  hazardous-vortex  encounters 
by  Heavy  and  Large  aircraft  following  3 nautical  miles  behind  Heavy 
aircraft.  The  probability  of  a hazardous-vortex  encounter  is  zero 
if  the  parameter  f,  the  fraction  of  the  roll-control  authority 
required  to  counter  the  rolling  motion  imposed  by  a vortex,  is 
about  0.5.  If  this  value  of  f is  acceptable  (higher  values  are 
accepted,  or  appear  to  be  accepted,  with  today's  safe  separation 
standards)  , then  the  VAS  may  be  used  to  reduce  Heavy/Heavy  and 
Heavy/Large  IFR  standards  to  3 nautical  miles  between  the  outer 
marker  and  touchdown  during  VAS-green  conditions. 

The  analysis  of  the  safety  of  3- naut ical -mi le  separations 
under  VAS-green  conditions  for  all  aircraft  required  a more  de- 
tailed model.  To  calculate  the  probability  of  a hazardous-vortex 
encounter  by  a Small  aircraft  behind  a Heavy  or  Large  aircraft, 
decay  alone  was  not  sufficient.  It  was  necessary  to  introduce  a 
cross-wind-aloft  model  which  contains  acknowledged  mathematical 
inconsistencies  but  yet  seemed  to  correlate  with  limited  data.  The 
analysis  also  required  comparing  conservative  estimates  of  hazard 
probabilities;  a possible  problem  with  this  method  was  indicated. 
Assuming  the  results  of  the  data  collection  (Volume  II)  and  the 
detailed  analysis  of  the  use  of  conservative  estimates  (Volume  III) 
validate  the  calculations  in  Section  5,  the  VAS  may  be  safely  used 
to  reduce  separations  between  the  outer  marker  and  touchdown  for 
all  aircraft  during  VAS-green  conditions. 


appendix  a 

INITIAL  VORTEX  STRENGTH 


The  strength  T of  the  vortices  shed  by  the  lead  aircraft 
parameterizes  the  degree  of  hazard  to  the  following  vortex- 
encountering  aircraft.  The  strength  or  circulation  of  a vortex 
may  be  used  to  determine  quantitatively  roll  moments  and  rates  an 
aircraft  would  experience  if  inadvertently  caught  in  a vortex. 

A Monostatic  Acoustic  Vortex-Sensing  System  (MAVSS)  (Refs.  9 
and  10) has  been  place  at  Chicago  O'Hare  on  the  approach  to  runways 
14R  and  32L  since  the  summer  of  1976.  The  function  of  the  MAVSS 
is  to  measure  a vertical  profile  of  the  vertical  component  of  the 
wind.  Since  the  vertical  component  of  the  ambient  wind  is  zero  at 
the  ground,  the  vortex  winds  are  measured  with  little  interference 
from  the  ambient  wind.  A single  MAVSS  antenna  measures  the  veloc- 
ity profile  of  a vortex  as  it  drifts  over  the  antenna;  an  array  of 
antennas  samples  the  vortex  at  different  times,  and  thereby, 
measures  vortex  decay. 

The  measured  vertical  velocity  distributions  are  used  to  cal- 
culate an  "effective  strength"  --the  strength  of  a line  vortex 
which  will  produce  the  same  torque  on  an  aircraft  as  the  measured 
velocities  will  produce.  In  other  words,  the  first  moment  of  the 
measured  vertical  velocity  distribution  is  defined  to  be  equal  to 
the  first  moment  of  a potentia1  or  line  vortex  (Ref. 10): 
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T'  is  the  "effective  strength,"  be  the  wingspan  of  the  vortex- 

encountering  aircraft,  and  V.,  , are  the  measured  vertical 

B Doppler 

velocities . 

At  Chicago,  ^p0ppjer  data  have  been  collected  for  1398R-707s, 

651 1 B- 72 7s , 758  B-737s,  495  B-747s,  744  DC-8s,  2797  DC-9s,  1073  DC- 
10s,  391  l. - 1 0 1 1 s , 223  B-707HS,  204  DC-8Hs,  356  small  propel  1 er - d riven 
aircraft,  and  68  Learjets.  The  reduced  data  consist  of  four  values 
of  T'--  P evaluated  for  semispans  (bg/2)  of  16.5  feet  (5  meters), 

33  feet  (10  meters),  66  feet  (20  meters),  and  100  feet  (30  meters). 

To  determine  rQ , the  initial  vortex  strength  required  in  the 
analytic  model,  the  maximum  values  of  P for  each  aircraft  type  and 
for  each  semispan  were  retrieved  from  the  MAVSS  data.  The  maximum 
values  were  used  as  the  stronger  the  vortex,  the  greater  the  hazard 
posed  by  that  vortex.  By  using  the  largest  observed  values  of  P, 
the  model  is  using  the  "worst  case"  in  the  sense  of  safety.  A 
linear  least-squares  fit  to  the  four  P values  (see  Fig.  17)  for 
each  aircraft  type  was  made  yielding  an  expression  of  the  form, 


Fo  = m be  + Fint  ’ 


where  m is  the  slope,  and  is  the  intercept  of  the  rQ  line  with 

the  P axis.  Note  that  there  are  many  data  points  not  shown 
below  those  plotted  in  Fig.  17;  for  example,  there  are  1397  data 
points  below  each  of  the  values  plotted  for  the  B-707--only  the 
largest  P value  was  plotted.  Thus,  the  initial  value  of  the 
vortex  strength  is  a function  of  the  wingspan  of  the  vortex-en- 
countering aircraft.  The  parameters  are  listed  in  Table  9 along 
with  the  correlation  coefficient  for  the  linear  regression. 
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TABLE  9.  PARAMETERS  FOR  T 

o 


Slope 

H 

Correlation 

Coefficient 

(ft/sec) 

(f t ^/sec) 

B - 70 7 

18 

920. 

6 

0 

996 

B-  727 

17 

95 

895. 

4 

0 

993 

B-  737 

11 

23 

642  . 

9 

0 

993 

B - 74  7 

19 

56 

1148. 

6 

0 

968 

DC-8 

13 

52 

1270. 

5 

0 

977 

DC-9 

12 

29 

837. 

0 

0 

998 

DC-10 

15 

46 

1010. 

3 

0 

969 

L-  1011 

16 

35 

958. 

5 

0 

977 

B-707H 

14 

88 

1007. 

4 

0 

966 

DC-8H 

16 

13 

982. 

2 

0 

990 

PA-28 

9 

54 

521. 

8 

0 

971 

Lear j et- 2 5 

5 

20 

715. 

2 

0 

999 

APPENDIX  B 


VORTEX  DESCENT 

The  purpose  herein  is  to  review  the  mechanisms  governing  vor- 
tex descent.  The  material  follows  closely  the  discussion  of  vortex 
transport  in  Reference  2. 


B.l  VORTEX  DESCENT  OUT  OF  GROUND  EFFECT 

For  vortex- transport  calculations,  the  Rankine  vortex  model 
is  usually  used: 


v(r)  = 


for  r < r , and 
c 


v(r)  = 


r 

2irr 


(32) 


(33) 


for  r > r . T is  the  strength  or  circulation  of  the  vortex,  and 
r^  is  the  radius  of  the  core.  The  basic  assumption  of  vortex 
transport  is  that  the  transport  velocity  of  each  vortex  is  the 
velocity  of  the  surrounding  velocity  field  evaluated  at  the  cen- 
troid of  each  vortex.  The  velocity  field  is  composed  of  the  am- 
bient wind  and  the  velocity  imposed  by  other  vortices.  Equation 
(33)  can  be  used  to  calculate  the  mutual  induction  of  vortices 
upon  each  other  if  rc  is  less  than  the  separation  distance  between 
the  vortices  (bQ)  . 

After  the  wing-vortex  sheet  has  rolled  up,  the  trailing  sys- 
tem consists  of  a pair  of  vortices  of  finite  rotational  core  area 
with  a core  radius  a relatively  small  fraction  of  the  vortex  span. 
If  the  vortex  pair  is  immersed  in  a "still,"  homogeneous  inviscid 
flow,  the  pair  is  convected  downward  at  a velocity  r/2TTbQ.  The 
classical  analysis  (Ref.  2)  shows  that  there  is  a closed  recircu- 
lating mass  of  air,  of  roughly  oval  proportions,  associated  with 


the  concentrated  vortex  pair,  and  that  the  cell  is  convected  down- 
ward at  a uniform  speed.  I-'low  exterior  to  the  cell  never  enters 
it.  Thus,  a long  vortex  pair,  which  may  be  regarded  as  substan- 
tially two-dimensional,  will  move  downward  in  an  unbounded  fluid 
with  constant  velocity  for  all  times. 

In  real  flows,  this  situation  does  not  persist  indefinitely, 
and  most  experiments  show  the  rate  of  descent  decreases,  and 
finally,  approaches  zero.  This  is  caused  by  diffusion  of  core 
vorticity  by  some  combination  of  laminar  and  turbulent  viscosity, 
and  it  will  occur  even  in  homogeneous  (unstratified)  flows.  Much 
effort  has  gone  into  explaining  and  quantifying  the  effect,  but 
the  subject  still  remains  controversial.  A rational  interpretation 
of  the  effect,  coupled  with  careful  observations,  which  greatly 
assists  explanation  of  the  effect,  has  been  put  forward  by  Maxworthy 
(Ref.  3 4)  . 

Maxworthy  conducted  experiments  with  vortex  rings  in  water, 
using  various  visualization  techniques  to  identify  where  the  flow 
went.  When  the  vorticity  was  relatively  well  distributed  in  the 
ring,  he  observed  that  the  outer  flow  was  entrained  into  the  back 
of  the  cell,  causing  an  increase  in  the  cell  volume.  At  the  same 
time,  a portion  of  the  cell  vorticity  was  shed  into  the  wake, 
removing  both  vorticity  and  momentum  from  the  cell.  The  combined 
effect  was  to  increase  the  cell  size  and  to  reduce  its  propagation 
velocity . 

The  mechanism  of  mass  entrainment  is  important  for  further 
development  of  vortex  transport.  Figure  18  is  a sketch  of  the 
vortex  flow  field,  in  coordinates  fixed  at  the  core  centers,  so 
that  the  outer  flow  is  represented  by  a uniform  and  unsteady  flow 
from  below.  The  cell  has  a well  defined  stagnation  point,  A,  and 
over  the  front  portion,  a well  defined  cell  boundary,  A-B.  Across 
this  boundary,  the  pressure  and  velocity  fields  of  inner  and  outer 
flows  are  continuous,  the  only  discontinuity  being  between  the 
inner  vortical  fluid  and  the  outer  irrotational  flow.  Because  of 
both  laminar  and  turbulent  effects,  the  vorticity  of  the  inner 
vortical  flow  is  transferred  to  the  outer  flow,  and  as  a 


consequence,  the  total  pressure  of  the  flow  is  reduced.  After 
passing  the  maxiinum-veloci  ty  point  near  B,  the  outer  flow,  con- 
tained approximately  by  the  stream  tube  DC,  is  unable  to  recover 
sufficient  velocity  to  rejoin  the  outer  flow  at  the  rear,  and 
remains  as  part  of  the  stationary  cell.  Thus,  the  cell  size  is 
increased.  At  the  same  time,  a neighboring  stream  tube  EF  acquires 
a smaller  amount  of  vorticity,  and  suffers  less  pressure  reduction, 
so  that  it  does  depart  from  the  cell  at  the  rear,  but  at  a lower 
than  free-stream  total  pressure.  This  portion  develops  into  a 
wake  behind  the  cell. 

The  same  process  causes  entrainment  of  the  outer  flow  into 
the  cell  and  a detrainment  (removal)  of  some  of  the  cell  vorticity 
and  momentum.  A further  process  occurs  on  the  centerline  of  the 
cell,  AX.  Here,  vorticity  is  annihilated  by  diffusion  from  the 
left-  and  right-hand  cells.  Thus,  three  vorticity  transfer  mech- 
anisms occur,  and  the  overall  effect  controls  the  cell  dynamics. 

Maxworthy  (Ref.  34)  showed  that  initially  the  vortex-shedding 
to  form  the  wake  was  extremely  weak  since  the  cell  vorticity  at 
the  boundary  was  quite  weak.  Thus,  although  the  cell  grew  in  size, 
it  did  not  lose  momentum,  and  the  impulse  was  conserved.  In  these 
circumstances,  the  main  vorticity  loss  occurred  along  the  center- 
line  and  was  small,  and  there  was  minimal  wake-momentum  loss. 

During  the  later  stages  in  growth,  when  more  vorticity  was  present 
near  the  boundary  between  the  inner  and  outer  flows,  the  wake 
developed.  Vorticity  and  momentum  were  shed  from  the  cell,  and 
thus,  the  momentum  in  the  cell  decreased  while  the  cell  size  in- 
creases. Both  of  these  effects  contributed  to  the  reduction  in 
descent  speed  and  the  final  complete  annihilation  of  the  cell 
momentum  for  the  two-dimensional  vortex  observed  by  Maxworthy.  It 
is  noted  that  for  a three-dimensional  vortex  pair  such  as  an  air- 
craft wake,  vortex  annihilation  caused  by  three-dimensional  effects 
(sinuous  instability  and  core-bursting),  almost  always  precedes 
the  two-dimensional  vortex  annihilation. 

It  must  be  noted  that  Maxworthy' s experiments  were  conducted 
with  vortex  rings  at  extremely  low  Reynolds  numbers  for  which  the 
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flow  was  certainly  laminar.  However,  f low-visual i zat ion  tests  with 
finite  wings  also  exhibited  a detrained  wake.  These  experiments 
were  also  performed  at  very  low  Reynolds  numbers.  It  is  possible 
that  during  the  later  stages  of  development  of  an  aircraf t- trai 1 ing 
vortex  system  that  similar  processes  of  mass  entrainment  and  momen- 
tum detrainment  occur.  For  laminar  transfer,  the  time  scales 
would  be  too  long  to  be  of  interest,  but  if  the  transfer  were 
assumed  to  be  turbulent,  it  might  be  possible  to  account  for  some 
of  the  observed  effects.  Thus,  it  appears  very  probable  that  the 
later  development  of  a vortex  pair  follows  qualitatively  the 
stages  described  by  Maxworthy  with  an  additional  initial  stage. 

The  three  stages  are  postulated  as  shown  in  Figure  19. 

Stage  I --  The  Inviscid  Cell:  Here,  the  vorticity  is  con- 
fined to  well  within  the  cell  boundary.  The  cell  boundary  is 
defined  as  the  streamline  between  the  flow  which  remains  with  the 
vortex  pair  and  that  which  remains  with  the  ambient  air.  On  the 
boundary  itself,  there  will  be  no  laminar  mass  or  momentum  trans- 
fer (since  there  is  no  distortion)  , and  turbulent  transfer  will 
have  no  net  effect  since  both  inner  and  outer  flow  have  the  same 
total  pressure.  In  these  circumstances,  the  core  size  is  less 
than  the  cell  size,  and  the  inviscid  cell  model  will  be  a good 
representation  of  the  dynamics.  The  time  rates  of  change  of  cell 
size  and  vortex  strength  are  zero;  the  descent  velocity  is  constant. 

Stage  II  --  The  Entraining  Cell:  As  the  core  vorticity 
diffuses  and  approaches  the  cell  boundary,  the  first  process  (of 
mass  entrainment)  occurs,  and  the  cell  grows;  F = 0.  The  descent 
velocity  reduces  slightly  from  the  inviscid  value  because  of 
greater  vortex  separation  resulting  from  cell  growth. 

Stage  III  --  The  Decaying  Cell:  During  the  later  stages, 
substantial  mass  entrainment  and  momentum  and  vorticity  shedding 
occur,  causing  a wake  to  develop  behind  the  vortex  pair.  Thus, 
various  catastrophic  instabilities  usually  develop  before  the 
complete  decay  has  occurred. 


FIGURE  19.  VORTEX  CELL  STAGES 


B . 2 VORTEX  BUOYANCY 


Vortex  buoyancy  is  the  aerostatic  force  imposed  on  the  vortex 
by  virtue  of  the  difference  in  density  between  the  air  contained 
within  the  vortex  and  the  surrounding  ambient  air.  There  are 
three  sources  of  this  density  difference:  The  first  is  a result 
of  the  static  underpressure  of  the  vortex.  The  second  is  the 
result  of  entrainment  of  hot  exhaust  gases  from  the  engines.  The 
third  is  the  result  of  descent  through  a nonadiabatic  atmosphere. 
The  first  two  effects  give  the  vortex  a positive  (upward)  buoyancy 
force.  The  third  may  give  a positive  or  negative  force  depending 
on  the  temperature-lapse  rate  of  the  surrounding  atmosphere.  The 
magnitude  of  these  effects  on  vortex  motion  is  small  compared  with 
mutual  induction. 

For  the  first  effect,  the  momentum  equation  for  cylindrical 
flow  gives 


(34) 


where  the  second  part  of  the  equation  is  valid  outside  the  core 
for  the  Rankine  vortex.  Assuming  adiabatic  conditions  in  the  core, 

~ k " k 

pp  = PooP<*>  = constant,  (35) 


where  k is  the  ratio  of  specific  heats,  and  the  subscript  00  refers 
to  ambient  conditions.  Equations  (34)  and  (35)  may  be  used  to 
find  dp/dr: 


k-1 


(36) 


Integrating  between  the  limits  of  infinity  and  r, 
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As  an  example  of  the  magnitude  of  the  effect,  for  the  flight 
conditions  for  a landing  B-727  and  an  assumed  core  radius  of  2 
feet,  the  density  is  0.00138  slugs/ft^  compared  with  an  ambient 
density  of  0.00235  slugs/ft^.  The  upward  buoyant  force  is 


F 


p (r))  rdr . 


(38) 


(In  practice,  the  integral  would  be  taken  from  the  core  radius, 
rc>  because  the  equation  is  based  upon  the  Rankine  vortex  assump- 
tion, and  is  only  valid  for  r _>  r^ . ) 

Engine-exhaust  extrainment  is  another  aircraft  variable  which 
can  affect  the  vortex-transport  process  through  variations  in  the 
vortex  buoyancy.  The  density  variation  caused  by  exhaust  entrain- 
ment is 


p = Poc  VT,  (39) 

where  T is  the  temperature  of  the  air  in  the  core  of  the  vortex, 
and  T^  is  the  ambient  air  temperature. 

The  third  cause  of  buoyancy  is  vortex  descent  through  a non- 
adiabatic  atmosphere.  The  interface  between  the  wake  fluid  and 
the  exterior  atmosphere  is  only  a dividing  streamline;  it  does  not 
support  shear  or  pressure.  In  the  presence  of  atmospheric  turbu- 
lence or  turbulence  in  the  wake,  the  streamline  is  perturbed,  and 
some  mixing  between  the  interior  wake  fluid  and  the  exterior  atmos- 
phere occurs.  When  the  mixing  is  small,  as  it  is  when  the  wake  is 
young  and  the  turbulence  is  low,  the  wake  fluid  retains  its  iden- 
tity and  its  physical  properties  as  it  moves  about  through  the 
atmosphere.  One  consequence  is  the  creation  of  aerostatic  forces 
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as  the  wake  moves  into  regions  of  varied  temperature  and  density. 
For  example,  an  upward  (or  buoyant)  force  is  developed  on  those 
segments  of  a wake  descending  into  a stably  stratified  atmosphere 
if  there  is  little  or  no  mixing  of  the  ambient  air  with  the  wake 
fluid.  Results  of  experiments  on  wake  buoyancy  are  presented  by 
Tombach  (Ref.  3 5);  lie  showed  that  a wake  descending  into  a stable 
atmosphere  acquired  buoyancy,  until  at  some  later  time,  turbulent 
mixing  between  the  wake  and  the  atmosphere  became  strong  enough  to 
erode  the  temperature  difference  between  them.  Increasing  ambient 
turbulence  shortened  the  period  during  which  buoyancy  was  acquired, 
and  as  a consequence  of  more  rapid  mixing,  resulted  in  more  rapid 
decay  of  both  buoyancy  and  vortex  descent. 

The  aerostatic  forces  are  a result  of  an  increase  in  the  tem- 
perature of  the  fluid  in  the  wake  oval , caused  by  adiabatic  com- 
pression as  the  oval  descends  into  a denser  atmosphere.  An  atmos- 
pheric temperature  stratification  other  than  adiabatic  (neutral) 
will  result  in  a temperature  difference,  and  hence,  a density 
difference  between  the  wake  and  the  atmosphere.  Buoyancy  will  thus 
be  created  for  wakes  descending  into  a stable  atmosphere  (the  most 
common  situation)  with  no  mixing.  Theoretical  models  with  no  mix- 
ing (Ref.  36)  indicated  that  the  buoyancy  so  acquired  accelerated 
wake  descent  and  decreased  vortex  spacing.  The  experimental 
observations  discussed  above  indicated,  however,  the  possibility 
of  a retarding  tendency  caused  by  buoyancy,  but  might  be  a con- 
sequence of  entrainment  rather  than  of  buoyancy. 

Overall,  the  effects  of  aerostatic  forces  on  vertical  wake 
motions  appear  to  be  of  smaller  order  than  the  dissipative  mechan- 
isms associated  with  turbulence,  which  could  overwhelm  the  buoyancy 
effect,  and  thus,  result  in  the  difficulties  experienced  in  prop- 
erly isolating  the  buoyancy  influence.  A comprehensive  discussion 
of  both  theoretical  and  experimental  observations  on  the  descent 
of  a wake  in  a stratified  fluid  is  given  in  Reference  36. 

The  predominant  effect  of  atmospheric  stability  appears  to  be 
the  indirect  one  associated  with  the  vertical  air  currents  result- 
ing from  atmospheric  mixing.  In  a stable  atmosphere,  this  mixing 

B - 9 


¥ 


is  suppressed,  resulting  in  reduced  vertical  air  motions  and 
reduced  effects  on  vertical  wake  motions.  In  unstable  conditions, 
vertical  atmospheric  activity  and  resulting  wake  motions  are  ampli- 
fied. Such  motions  may  be  either  upward  or  downward,  thus  unpre- 
dictable wake  ascent  or  descent  can  result  under  such  unstable 
atmospheric  conditions. 

The  effects  of  atmospheric  stratification  on  initial  wake 
descent  rates  are  shown  schematically  in  Figure  20.  For  a stable 
atmosphere,  the  wake  descends  initially  at  a speed  which  is  con- 
sistent with  the  inviscid  analytical  model.  Random  influences 
become  more  evident  for  less  stable  conditions.  The  figure  applies 
only  to  vertical  wake  motion  during  the  first  few  moments  after  the 
wake  has  been  fully  formed.  Subsequent  vertical  wake  motions  are 
influenced  by  buoyancy,  turbulence,  and  the  continued  random  action 
of  vertical  air  motions  (which  are  accentuated  for  less  stable  con- 
ditions) . Not  all  of  these  factors  can  be  quantified  at  present. 

Vertical  air  motions  often  become  quite  pronounced  under  the 
influence  of  thermal  activity  near  the  ground.  Extreme  vertical 
convolutions  of  wakes  can  occur  because  of  convection  by  these  air 
motions.  Under  these  conditions,  the  vortex  pair  is  stretched  and 
distorted  into  a highly  nonlinear  configuration  where  mutual  and 
self-induced  vortex  velocities  act  to  amplify  the  process.  Figure 
21  shows  such  a wake  generated  near  the  ground  under  conditions  of 
high  atmospheric  thermal  activity.  The  impossibility  of  dealing 
with  such  a resultant  wake  structure  on  anything  but  a statistical 
basis  is  clearly  evident.  However,  it  remains  to  be  shown  that 
such  a convoluted  wake  may  pose  a hazard  in  operational  conditions. 

All  wake  motions  near  the  ground  do  not  exhibit  such  extreme 
behavior.  Under  stable  atmospheric  conditions  and  reduced  thermal 
activity,  the  wake  and  the  vortex  pair  undergo  more  orderly  motions 
which  are  fairly  well  understood  and  can  be  approximated  analyti- 
cally. luch  conditions  are  also  those  of  greatest  operational 
interest  because  these  same  factors  are  conducive  to  wake  persis- 
tence. Wakes  generated  or  moving  into  ground  effect  are  subjected 
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to  the  influences  of  the  induction  velocities  generated  by  the 
underground  or  image  vortices.  The  vertical  motion  of  the  wake 
gradually  and  predictably  slows  and  eventually  stops  at  an  alti- 
tude equal  to  about  one-half  the  initial  vortex  spacing.  The 
vortices  then  move  apart. 

B . 3 DESCENT  DATA 

Tombach  (Ref.  3 5)  studied  the  descent  of  vortices  shed  by  a 
light,  twin-engine  aircraft  (AeroCommander  560F)  . Cameras  recorded 
the  motion  and  decay  of  smoke-marked  vortices. 

Descent  trajectories  for  the  AeroCommander  vortices  are 
shown  in  Figure  22  for  unstable  and  neutral  stratifications,  and 
in  Figure  23  for  stable  isothermal  and  inversion  conditions.  The 
solid  line  in  the  figures  is  the  theoretical  descent  rate  based 
on  elliptical  wing  loading.  The  broken  lines  indicate  other  de- 
scent rates  for  comparison. 

In  Figure  22,  the  wakes  actually  rise  in  the  two  unstably 
stratified  cases  shown,  probably  due  to  being  carried  upward  by 
the  considerable  vertical  currents  which  accompany  instability. 

The  high  turbulence  which  naturally  occurs  in  such  an  unstable 
atmosphere  results  in  very  brief  lives  for  these  wakes.  The  wakes 
in  a neutral  atmosphere  show  a fairly  rapid  descent,  with  initial 
speeds  well  exceeding  the  theoretical  descent  rate.  After  20  to 
30  seconds,  the  descent  has  slowed,  but  the  speeds  still  exceed 
the  theoretical  rate. 

For  the  more  stable  atmospheres  (Fig.  23),  the  range  of 
initial  descent  speeds  bracket  the  theoretical  descent  rate  and 
are  within  about  25  percent  of  it.  As  before,  there  is  often  slow- 
ing down  of  the  speed  of  descent  after  30  to  40  seconds,  with 
descent  speeds  at  50  seconds  typically  about  1/2  to  3/4  of  their 
initial  values. 

Wakes  often  exhibit  considerable  banking  or  tilting  of  the 
plane  containing  the  vortex  pair.  Figure  24  displays  trajectories 
for  wakes  which  have  banked  in  excess  of  25  degrees.  Surprisingly, 
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FIGURE  24.  DESCENT  TRAJECTORIES  FOR  BANKED  WAKES.  The  theoretical 
descent  rate  for  unbanked  wakes  is  indicated  by  the  solid  line. 


banking  does  not  change  the  descent  speeds  greatly;  the  trajecto 
ries  do  not  look  much  different  than  those  for  unbanked  wakes. 

The  data  show  the  effects  of  atmospheric  stability  on  the 
descent  of  vortices.  Slowing  of  descent  appears  to  be  caused  by 
the  turbulent  erosion  of  the  wake  strength  (with  an  attendant 
decrease  in  the  potential  hazard  if  the  wake  were  encountered  by 
an  aircraft) . The  data  are  consistent  with  the  three-stage  char 
acterization  of  the  life  of  a wake  discussed  in  Section  B.l. 
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APPENDIX  C 


CROSS-WIND  PROBABILITY  DENSITY 

C.l  GENERAL  DERIVATION  OF  CROSS-WIND  PROBABILITY  DENSITY  USING 
WIND -ROSE  DATA 

Wind-rose  data  are  ordinarily  presented  as  a polar-coordinate 
p.d.f.  (or  more  usually  as  a discrete  probability  distribution) 
for  the  joint  probability  of  a wind  magnitude  (r)  and  "from"  direc- 
tion  (9).  Denote  this  wind-rose  as  p (r,0). 

To  determine  the  p.d.f.  for  the  cross  wind  on  a given  runway, 
consider  the  following  transformation  of  variables,  where  a is  the 
runway  heading,  and  (x,y)  defines  a cartesian  coordinate  system 
with  the  positive  x-axis  aligned  with  the  runway  heading  (Fig.  25). 
The  transformations  from  the  wind-rose  to  the  runway-coordinate 
system  are: 


head  wind:  w = r cos(a  - 0) , 

x 

cross  wind:  wy  = r sin(a  - 0), 


and  the  inverse  transformations  are: 


2 2 2 

magnitude:  r = w + w , 

x y 

direction:  0 = a - tan  * (w  /w  ) 

y x 


The  standard  formula  for  transforming  the  joint  p.d.f.  is 


p (w  ,w  ) = 
*wv  x y 


xy  1 


p!  (CZ 


T 

w +w 
X y 


a - tan 
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(wy/wx)),  (40) 


where  |J 

1 xy 

formation : 


is  the  absolute  value  of  the  Jacobian  of  the  trans- 
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Thus , 


(42) 


is  the  representation  of  the  wind  p.d.f.  in  cartesian  coordinates 
aligned  with  the  runway  heading  a. 

Since  the  concern  in  vortex  analysis  is  with  the  cross-wind 

magnitude,  define  wc  = | w^. | . The  marginal  p.d.f.  for  wc  is  found 

by  integrating  out  w at  w = +_  w : 

x y c 


MV  = / 'PwM-V  + Pw(V-wcmV  (43) 

- 00 


At  this  point,  it  is  useful  to  make  a simplifying  assumption  that 
the  wind  distribution  is  symmetric  about  the  x-axis;  i.e., 


Then 
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Pc(wc)  = 2 / Pw(wx,wc)  dwx> 
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dw 


(44) 


^wx  + WC 


giving  the  cross-wind  p.d.f.  from  the  wind  rose. 

Wind-rose  data  are  a statistical  representation  of  the  proba- 
bility of  a particular  wind  vector  being  randomly  encountered  on 
(or  above)  the  airport.  The  transformation  to  cross  wind  on  a 
given  runway  is  strictly  speaking  only  valid  for  a randomly  se- 
lected runway.  Selecting  a particular  runway  for  landing  is  an 
act  implying  prior  knowledge  of  the  wind  magnitude  and  direction, 
the  usual  strategy  being  to  pick  a runway  heading  (a)  to  minimize 
tail  and  cross  winds.  Under  such  conditions,  limits  on  the  allow- 
able cross  and  tail  winds  should  be  imposed  when  defining  the 
cross-wind  p.d.f.,  given  that  runway  a has  been  selected.  For 
example,  let  be  the  maximum  allowable  tail  wind,  and  Lc  the 
maximum  cross  wind,  then  an  improved  estimate  of  the  cross-wind 
p.d.f.  when  runway  a is  in  use  would  be 


00 

/p  (w  ,w  ) dw 
L x’  cJ  x 


J / l w 
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/ C I p (w  ,w  ) dw  dw 
JQ  J_L  ‘w'-  x’  cJ  x c 


(45) 


This  will  be  further  modified  by  runway  usage  and  noise  considera- 
tions . 
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C . 2 CROSS-WIND  PROBABILITY  DENSITY  DERIVED  FROM  TOWER  DATA 


Crawford  and  Hudson  (Ref.  31)  measured  winds  using  acrovanes 
mounted  on  a 1500-foot  tower  located  near  Oklahoma  City.  Reference 
31  presents  annual  relative  frequency  distributions  of  wind  speed 
in  2-knot  intervals  for  various  levels  on  the  tower.  The  concepts 
of  Section  C.l  will  now  be  applied  to  the  wind  data  of  Crawford 
and  Hudson  to  obtain  a transformation  to  cross  wind. 

A number  of  assumptions  are  required  to  derive  the  cross- 
wind  p . d . f . : 

A1 . The  magnitude  of  the  wind  vector  is  as  described  in 
Reference  27. 

A2.  The  wind  direction  is  equally  likely  in  any  direction 
(i . e. , random) . 

A3.  Wind  magnitude  and  direction  are  independent  random 
variables . 

A fourth  assumption  is  implied  by  this  approach,  as  discussed 
at  the  end  of  Section  C.l; 

A4 . The  selection  of  runway  direction  is  independent  of  the 
wind  behavior  above  the  surface. 

Crawford  and  Hudson  state  that  the  wind  magnitude  measured 
about  1,000  feet  above  ground  could  be  well  approximated  by  a form 
of  Rayleigh  distribution 
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r > 0 


(46) 


where  a is  related  to  the  mean  observed  wind  speed  y by 


o 


y . 


(47) 


By  A1  above,  this  Rayleigh  distribution  can  be  taken  to  represent 
the  marginal  p.d.f.  for  the  speed  of  winds  aloft  at  the  airport. 


By  A2,  the  marginal  p.d.f.  for  wind  direction  can  be  written 

P0(0)  = ^ , (48j 

for  0 _<  8 _<  2 tt  . Combining  them  with  A3,  the  assumed  wind  rose  for 
winds  aloft  becomes 


2 , 2 
-r  /2o 


(r,6)  = 


2iTa 


(49) 


for  r ^ 0 and  0 _<  0 _<  2 tt  . A4  states  that  this  wind  rose  holds 
without  regard  to  the  runway  in  use;  without  loss  of  general ity, 
the  runway  direction  may  be  disregarded,  and  a may  be  set  equal 
to  zero.  The  cross-wind  p.d.f.  is,  then, 


J (W  ) = -±y  f 

C C T TO2 


, 2 2,  2 

(wv  + w )/2o 


dw 


2/t  2 

2 'wc  /2° 


(50) 


with  wr  _>  0.  Note  that  this  is  the  p.d.f.  for  the  magnitude  of  a 

Normal  random  variable  (i.e.,  w is  distributed  Normally  with  zero 

2 c 

mean  and  variance  a ),  which  is  to  be  expected,  since  the  Rayleigh 
distribution  is  the  p.d.f.  for  tiie  vector  sum  of  two  identically 
distributed  Normal  random  variables.  Falls  and  Brown  (Ref.  35) 
substantiated  that  wind  components  may  be  considered  bivariate 
Normally  distributed. 

The  cross-wind  p.d.f.  may  be  rewritten  in  terms  of  the  mean 
wind  speed  aloft  (u) 


Pc(wc) 


2 2 

1 -ttw  V4p 

— e c 

y 


w > 0 
c — 


(51) 
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and  the  mean  cross-wind  magnitude  is  related  to  the  mean  total 
wind  by 

2 

E(wc)  = £ V . 

The  observations  reported  by  Reference  31  give  p = 18.6  knots  at 
1166  feet  above  ground  level,  with  17.2  and  18.7  knots  at  the  next 
lower  (873-foot)  and  higher  (1459-foot)  levels,  respectively. 

Using  the  18.6-knot  value,  the  derived  p.d.f.  for  cross-wind  mag- 
nitude is 

Pc(w  ) = 0.0538  e-°-00227w  2 
c c c (52) 


! 
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with  a mean  of  11.8  knots. 

If  the  winds  of  interest  are  correlated,  even  weakly,  with 
runway  selection,  then  factors  come  into  play  that  could  invali- 
date A2  to  A4  above,  and  the  more  specific  analysis  must  be  applied. 
The  effect  of  selecting  a specific  runway  for  approaches  into  the 
wind  will,  in  general,  further  limit  the  variance  of  the  cross 
wind,  making  low  cross-wind  speeds  even  more  likely.  Acknowledging 
this  effect  in  the  analysis  could  have  distinct  effects  on  the 
safety  calculations  for  vortices  of  non-zero  potential  hazard 
since  such  calculations  rely  in  the  main  on  cross-wind  transport 
effects  for  reductions  in  encounter  probabilities. 

C . 3 CROSS-WIND  PROBABILITY  DENSITY  FOR  A VAS-GREEN  LIGHT 

The  procedure  is  much  less  clear  when  attempting  to  infer  a 
cross-wind  p.d.f.  aloft  given  a "red/green  light"  observation  of 
surface  winds.  If  there  were  a simple  one-to-one  mapping  between 
the  surface  wind  vector  and  the  wind  vector  at  altitude,  then  an 
appropriate  determination  of  the  wind  rose  aloft  conditional  on  a 
surface  "green  light"  could  be  made.  However,  it  is  not  reasonable 
to  expect  that  such  a mapping  can  be  established,  or  that  any  sim- 
ple one-to-one  correspondence  should  exist.  In  fact,  it  is 
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reasonable  to  expect  that  a statistical  rattier  than  deterministic 
relationship  will  exist--that  is,  any  particular  wind  vector  ob- 
served on  the  ground  may  relate  probabilistically  to  a large  set  of 
possible  wind  vectors  aloft  (and  vice  versa).  This  makes  the  aloft 
interpretation  of  a "green/red  light"  surface  observation  much  more 
cumbersome  and  its  description  dependent  on  observed  statistical 
data--the  measured  wind  rose  aloft  given  a "green  light"  at  the 
surface.  The  latter  might  indeed  bear  no  relationship  to  surface 
winds,  in  which  case  the  observation  of  a "green  light”  would 
convey  no  information  of  use  in  establishing  the  p.d.f.  for  cross 
winds  aloft. 

In  the  absence  of  specific  meteorological  rules,  a simple 
assumption  is  postulated: 

A5.  If  the  surface  wind  is  observed  outside  the  VAS  ellipse, 
then  the  wind  vector  aloft  will  also  alway3  lie  outside  of  such  an 
ell ipse . 

This  assumption  has  several  drawbacks  and  further  assumptions 
implicit  in  its  use:  Steering  of  the  wind  vector  may  cause  a vec- 
tor nominally  outside  the  ellipse  near  the  MM  to  be  inside  the 
ellipse  when  applied  aloft.  A given  surface-wind  observation  may 
yield  a probability  distribution  aloft  that  is  radically  different 
from  the  normal  wind  rose;  when  an  area  (such  as  the  "green  light" 
region)  of  the  surface-wind  distribution  is  considered,  the  re- 
sultant conditional  wind-rose  p.d.f.  may  not  be  a simple  scaling 
up  of  the  ground-wind  rose,  and  most  certainly  will  not  look  much 
like  a scaled  surface-wind  rose.  Thus,  caution  is  needed  in  the 
use  of  A5. 

A further  assumption  is  implicit  in  the  application  of  AS: 

A6 . If  a "green  light"  occurs,  winds  aloft  maintain  their 
normal  (relative)  probability  relationship  that  was  expressed  by 
the  unconditional  aloft  wind  rose. 

This  results  in  a rather  unique  mapping  between  the  surface  and 
aloft--not  quite  deterministic,  but  extremely  rigid.  It  implies 
sufficient  correlation  between  the  two  wind  vectors  to  be  at 
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direct  odds  with  A2-A4.  This  remains  an  unanswered  contradiction 
to  be  resolved  by  the  results  of  the  data  collection  (Volume  II  of 
this  report).  With  A5  and  A6 , the  determination  of  the  conditional 
p.d.f.  for  winds  aloft  given  a "green  light"  at  the  surface  is 
straightforward . 

The  VAS  ellipse  (E)  is  oriented  with  the  runway-coordinate 

system;  the  semi-major  axis  a^  of  the  ellipse  is  aligned  with  w^, 

and  the  semi-minor  axis  bF  with  w . Define  the  conditional  cross- 

r y 

wind-aloft  p.d.f.  given  a surface  "green  light  on  runway  " obser- 
vation to  be  p (w  /GLa) . This  is  found  (see  Figure  26)  by  inte- 
grating the  unconditional  aloft  p.d.f.,  Pw(wx>wy)>  along  the  two 
lines  Wy  = + wc , and  then  normalizing  by  the  prior  probability 
that  the  wind  vector  (w)  is  outside  the  ellipse,  p(w^E) . (Note 
that  this  is  not  the  a priori  probability  of  a "green  light,"  which 
is  the  same  calculation  applied  to  the  surface  wind  rose.  In  gen- 
eral, when  applied  at  altitude,  the  value  obtained  will  be  larger 
than  the  VAS  effectiveness  at  the  surface  because  of  the  higher 
wind  speeds  expected  at  altitude,  some  of  which  in  fact  may  occur 
with  "red  light"  conditions  at  ground  level.)  Thus, 


Pc  Oc/GLa) 


/[ p (w  ,w  ) + p (w  , -w  ) ] 

L ^W1"  X ’ c'  *WV  X ’ CJ 


dw 


p(w^E) 


00 

/[p  (w  ,w  ) + p (w  , -W  ) 1 

llWV  X ’ CJ  *WV  X CJ  ‘ 


dw 

X 

p(w^E) 


(S3) 


where 


0 i V bE  ’ 

br  < w < L 
E—  c — c 
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North 


aza>  Lines  of  constant  cross  wind 


"Green  Light"  Region 
(normalizing  factor) 


L I CURL  26.  R LG  IONS  Of  INTEGRATION  FOR 
"VAS  GREI.N  L I (.H  I ON  RUNWAY  •" 


As  before,  simplifications  can  be  made  regarding  the  symmetry  of 
cross  wind,  and,  when  data  show  no  prevailing  winds  from  any  direc- 
tion, then  head-tail-wind  symmetry  might  be  acceptable  also,  yield- 
ing 


an 

pR(\/w  2+  wr2,a-tan‘1(w  /w  )) 
w » x c 

dw 

x , 

(54) 

pw(wc/GLa)  = 4 

•V  / 2 2 

Vwx  + Wc 

p (wf^E) 

with 

OO  oo 

p(wf!E)  = 4 J J 

Pw(Vwx2+  wc^  a-tan’^W3 
Vwx2  + Wc2 

dw  dw 
x c 

• (55) 

Since  symmetry  in  all  directions  has  been  in  effect  assumed,  a can 
be  set  to  zero  without  any  further  loss  of  generality. 

The  assumptions  in  the  application  of  Pu,(w  /GLci)  are  A5  and 
A6  plus  any  conditions  that  guarantee  that 

A7.  Tiie  unconditional  aloft  wind  rose  is  symmetrical  with 
respect  to  the  selection  of  runway  direction. 

Note  that  condition  A7  (symmetry)  is  satisfied  by  A2  and  A3 
(independent,  equally  likely  wind  direction),  but  that  the  same 
claim  of  lack  of  surface/aloft  correlation  that  was  stated  as  A4 
still  remains  at  odds  with  the  basic  physical  premise  of  using  a 
surface  observation  to  provide  information  on  wind  behavior  aloft. 

The  assumed  Rayleigh  wind  rose,  converted  to  cartesian  runway 
coordinates,  is 


p (w  ,w  ) = 
^w^  x’  y 


Pw  CV*x2  * WvZ>  tt-tan-^Wy/Wx)) 


/w  + W' 

x y 
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! -(wx2  ♦ wy2)/2o2 

2 e 
2tio^ 


Thus,  by  substitution. 


Pw(wc/GL«)  = 


-(wx2  ♦ wc2)/2o2 


_0O 

XX 


-(wv  Z + w )/2o" 


[e/o)e 


X c 


2,,  2 

-wc  / 2o 


dwxdwc 


_CO  L , n 

r -w  /2a 

Jo  Q(e/o)e  dw£ 


where 


00  1 

Q(e/o)  = fz/g  ^ 


2/-, 
-v  /2 


=vr* 


(y  is  the  mean  local  wind  aloft)  and 


E L 2 ,2 

VbE  ' “c 


0 < w < b„ 
— c E 


w > b_ 
c — E 


for  the  conditional  p.d.f.  for  cross-wind  magnitude  aloft  given  as 
a VAS  "green  light  on  runway  a"  and  A1  to  A6  are  valid 
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(notwithstanding  the  unanswered  contradictions  among  the  assump- 
tions). Given  the  parameters  jj=  18.6  knots,  3j.  = 12.5  knots, 
bp  = 5.5  knots,  and  cx  = any  runway  direction,  the  equation  can  be 
evaluated  numerically,  and  is  shown  in  Figure  13  along  with  the 
no-VAS  condition.  The  VAS  "green  light"  observation,  as  expected, 
implies  a shift  in  the  p.d.f.  from  lower  to  higher  cross-wind 
values . 

A more  complete  discussion  of  the  material  in  this  appendix 
can  be  found  in  Reference  38.  In  particular,  the  assumptions  and 
the  implications  of  the  assumptions  are  discussed  in  detail. 

C.-1  DATA  COLLECTION 

When  the  surface  wind  is  outside  the  VAS  ellipse,  does  the 
wind  vector  aloft  also  lie  outside  the  ellipse?  This  is  an  assump- 
tion (A5) , which  was  invoked  in  the  absence  of  any  specific  mete- 
orological rules  relating  winds  at  the  surface  to  winds  along  the 
ILS  out  to  the  OM.  To  test  the  assumption,  wind  data  extracted 
from  Aircraft  Integrated  Data  System  (AIDS)  recordings  were  used. 

The  AIDS  package  consists  of  an  onboard  computer  which  records 
Inertial  Navigation  System  data  along  with  aircraft-control  inputs 
(throttle,  flaps,  etc.).  The  data  came  from  62  KLM  B-747  landings 
at  Chicago  O’Hare  International  Airport.  The  wind  was  determined 
in  approximately  25-foot  intervals  between  50  and  1300  feet  above 
ground  level.  The  flights  occurred  between  1 July  and  23  October 
1976. 

Whenever  the  wind  exceeded  10  knots  at  the  surface,  little 
turning  or  change  in  wind  direction  was  noted  among  the  various 
altitudes  (usually  less  than  20  degrees).  Between  6 and  10  knots, 
the  wind  direction  varied  up  to  40  degrees.  Below  6 knots,  large 
variations  in  wind  direction  were  noted  (sometimes  exceeding  100 
degrees).  Fifty  of  the  62  approaches  occurred  during  "VAS-green” 
conditions.  For  the  50  approaches,  the  winds  aloft  were  in  agree- 
ment with  A5  for  48  approaches  (96  percent).  One  of  the  two 
approaches  not  conforming  to  A5  yielded  a discontinuous,  and  hence, 
suspect  wind  profile,  whereas  the  other  case  may  have  occurred  in 


j 
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"VAS-red"  conditions  as  the  surface  wind  was  on  the  edge  of  the 
VAS  ellipse  (a  5.5-knot  cross  wind).  However,  it  does  appear  that 
AS  is  warranted. 

It  is  acknowledged  that  the  AIDS  data  are  restrictive  as  the 
B-747  landings  occurred  at  about  the  same  time  each  day  (approxi- 
mately 1600  Local  Standard  Time).  Until  the  completion  of  the 
data  collection  discussed  in  Appendix  F (the  analysis  will  be  con- 
tained in  Volume  II  of  this  report),  the  AIDS  data  are  the  only 
relevant  data  available;  they  don't  contradict  A5  or  the  cross- 
wind  aloft  model . 
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LISTING  OF  COMPUTER  PROGRAM 

The  encounter-probability  model  described  in  Section  4 was 
coded  in  FORTRAN  for  running  on  the  TSC  DECsystem- 1 0 . Figure  27 
is  a listing  of  the  routine.  Many  comments  have  been  inserted  to 
permit  identifying  the  corresponding  equations  of  the  probability 
model.  Lower  case  letters  cannot  be  used  in  programming  the  DEC 
system-10.  Usually  the  notation  used  in  the  program  is  a straight- 
forward interpretation  of  the  parameter  in  the  text  (e.g.,  F for 
f,  SIGD  for  Op,  T or  TT  for  t,  etc.).  Examples  where  the  program 
notation  may  be  confusingly  inconsistent  with  the  text  are  listed 
below : 


TEXT 

b 


D 

U 


Pc(wc) 


PROGRAM 

BF 

B 

BG 

DHAZO 

EL 

PCW 

RR  or  R 

SIGWD 


D-l 


t t 


C...  F.NCCuNTtR  PR'.bAulL  1 T y ANALYSIS 
C...  0(t  T / T SC  u.klNSTrN  CuOl  521 

0 OMMON/P  C.v  INU/PCM  5 1,2) 

DIMENSION  KlCCPCc)  , 1 VAS ( 2 I 

L'lMtNS  lot)  PJESUM(o)  , III  ol  , PJE  l 01  ,RR  16)  ,UOIo>  ,P0t  »(  61 
OIMFNSKN  S1GCC16) 

C...  Nil  VAS  , VAS  CRuSSMNO 

DATA  I VAS/  • NO*  , • ».  I TH  •/ 

DATA  NLUUP/1,1/ 

UPO  M JN  1 T = 5,  DlV  I C E=  *LPT  • , ACCESS^  • SlUOUT  • .MOufc  = • ASCII  • . 

1 F I L E = ' £ NC  I k . * I 
CALL  PkSET 

10  C L A T 1 ACF 

C • • • £ IS  THt  FKAC  T ION  OF  kLLL“Cl.\TkiiL  AuTHUklTY 

F=0. 379 

0*** 

or  loco  jjj=  l,  i ) 
j = jj  j 

1 F ( J .EJ.l  i ) CD  TO  AOO 

or  loco  jFj=i(iL 

JF= JFJ 

CALL  C AT FIL(J,IAl,B,V, SLOPE, BINT, VC, SI CMD) 

CALL  CATF  1L1JF  ,1ACF  *bF,U,SLOPOF,BINTF  ,VL)F,S1GMlF  1 
C...  T 1 PE  INCPfcKcNT  C SEC/ MM  » 

DT=fcC7E./V 

C...  ulSTANCt  (FT)  F°CM  RUNwAY  THPtSHOLC 
tL=3MCCC.O 
JL=  ROCO.O 

C...  GAMMA  FIT 

GAM=(SLQPl*(6F/3.28J8AI*61NTI*(3.28J8mI**2 
k F I T E ( 5 , AOOO ) F, I AC, V.O.GAM, VU.S IuMJ.SLOPl , BINT 
ACOO  FORMAT  ( IH1  , • ENCCcMEk  ANALYSIS*  ,3X,*F=*,F5.A/ 

* • CfcAEFAI ING  A/C  •/ IX, AS,  1 X,  *V  = * ,F fc.  1,  1*, *B  = • ,F6.  1 ,1  A, 'UAMMAC  =• 

* ,F6.1,1X,'VJ  -•  , F 6.2/ 5X  , 1 S IoMJ  =•  .Fu.z.lX,  'SLCPt  *• 

* ,Ffc.  3,  IX,  • INTERCEPT  *',F7.3/» 

WP  i T t ( 5 ,SOOO  1 lACF.U.KF 

500C  FORMAT  < lx,  • FOLLUxcR  A/C  • / IX ,rt 5, 2X , *0  * * ,F  6.  1 , 1 X , ■ BF  = *,F6.l) 
C...  PERFORM  ANALYSIS  FCk  FOLLL’kFR  OISIANCFS  UP  in 

C...  6 NM  AT  1-NPI  lNChFMtNTS  BEHIND  GENEkAT  CR 

C.  .. 

C...  CALCULATE  MINIMUM  DISTANCE  FCk  0 PROBABILITY 
UHAZOMI  GAM  *8  )/  ( BF  *U ) )*(  20.  52o* 9 . 5K/6C7&.  ) /F 
lF(JF.ut.ll)OHAZO=C.75*OPAZO 
lDHu,  = C 

000=97.9*8/6076.0 
IF  (CFA20.LT.0Uul  IDHC=  1 
IF  (CFAZO.LT. DCU)  CHAZ0  = CL)D 
„RITE(5«5CC1)  JHAZO.IOhU 

50  C l FORMAT  ( / * MINIMUM  C1STANCE  FOR  0 PRCBAiil  L I TY  ■ , F 7.2  , • NM.'.IT) 

T = 2. *0  T 
DC  250  1 = 3,6 
T=  T*D  T 
TT I I 1 = T 

C...  CALCLLAT ICN  OF  E NCCUNT  ER  CROSS  SEOTIUN 
0=  ( 0 . 2 * V *T  1/3 
GAMMA  = GAM 

C..  GAMMA  ATTENUATION  AFTEK  BREAKPOINT 

IF  I C.GE  .9. SB  I uAMMA  = 9.58*GAMMA/0 
C...  ENCOUNTER  CROSS  SECTION  (Rl 
R P ( I ) = ( 2 .653* G AMM A 1 / ( U*F I 
IFIJF.GE.il  > RR(  I I =0.  75*RM  I ) 


FIGURE  27.  LISTING  OF  COMPUTER  PROGRAM  OF 
PROBABILITY  MODEL 


THIS  PAGF  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURtUSHfiO  XO  DLQ 


C...  VCRTEX  DESCENT  PARAMETERS 
DOT  I ) = VURT 
SIGUD1  I)=S1GMU»T 
2 50  CON  TINGE 

E L = E L ♦ DL 

IE  lEL.ub. 50000.0)  GO  TO  1000 
SlGY=C. 011208* EL-5. A20E3 
SIGZ=C.003d6b*ELY9.d0A895 
wKITE(5,50C0)  El.SIGY.SIGZ 

9000  FCPMATl//'  L = 1 , E8 . 0 , ' F EET ■ , 5X , • S I GM Y= • , E6 . 0 . 5X , • S I GMZ = • , E fc. C t 
DO  1 500  KL  1 = 1 ,2 
DO  255  [=5,o 
255  POE  SEP  (I  ) =0 

KKL1=KL(jUPIKLI ) 

LL=KKl 1-1 

RR1 TE ( 5 ,0000)  lVAS(KLi),LL,([,I=3,0),(TT(l)«l=3,0) 

0000  FC.EMAT  I/A5,*  VAS'/'  PDESUM',13,'  - 50  K T S ■ /A  I 5 X , I 1 , • NM',5X)Z 
1 Al  2X,F  0.  1 ,'  SECS'  ) ) 

RrITE(5,6000)  ( RR 1 1 ) , 1 = 3,0) 

0600  rCKMAT  U13X.E6.1,'  FT  E T • ) ) 

C...  PIRFLRM  ANALYSIS  FOR  EACH  CKUSS-rINu  SPEED  AND  ACCUMULA  TE 

C...  PRCBAE 1L 1 TY  AS  A FUNCTION  CE  CRCsS-hIND 

DC  300  K=  KKL 1 , 5 1 
CC  A 1=3,0 
P = R R I 1 ) 

PCE  ( I ) = 0. 

PDER  ( 1 ) = 0. 

HAZEi  = R-BF/2. 

IEIHAZH.LE.O.)  jC  IP  A 
K0  = 0 

bG=3.C»(3.1A159l»E/B.O 
IE  (R.LT.3U)  Kr  = 1 
IF  (F.LT.du)  R = 0./07*R 
D=OU  1 I ) 

T=TT  I 1 ) 

SIGD  = SICDD(  I 1 
C...  WIND  RUN 

R=  E L C A T (K-l ) *1 .60b*T 
C...  DEVIATION  IN  RIND  RUN 
SICRD=15.0*T*l.O8b 
C...  HC-KlZCNTAL  deviation 

S I G E*  = SCCT  1 2 . 0 * S I OY*S  IGY  ♦ 5 I C„0*S  1 G = u ) 

C...  VERTICAL  DEVIATION 

SIGV=SCRT(2.0*SICZ*SluZ*oIuU*SIuD) 

c...  calcclatf  horizontal  e no c uni r r probability  irehi 

CALL  ERFUNU.HAZH.SIur.PEM) 

C...  LALCLLAT  E VERTICAL  ENCLuNTER  PROBABILITY  (PFV) 

HAZ  V = SCRT I hAZH* (R*bF /2. ) ) 

CALL  EPEUN(U,HAZV,SIGV,PEV) 

C...  PROBABILITY  OF  A hAZARDCLS  tNCCUNTER  (PCE) 

PCE ( l)  = PEH*PEV 

c...  rhen  f less  than  Eo  f,ut  uRlater  than  bf/z,  have  tc  cunsicer 

C...  EACH  VUKTEX  SEPARATELY 

IE  (KC.EC.l)  PDF ( I ) =2 .J*pnc ( I I 

C...  PROP  AT 1L 1 TY  OF  A HAZARDOUS  tNCucNTEK  ulVEN  THE  CPuSS-RIND 
PUT  R ( I ) = PDF  l I I *PC» ( R , K L I ) 

c...  probability  of  a hazaruulS  encllnier  accumulated  up 
C...  TC  A CRDSS-UNl  if  GIVFN  -AuNITUOE 

PDESUM( 1 l=PJESUM( 1 TRHJcrI I ) 

a cr:N  tinge 

300  CC  NT  I NGT 


FIGURE  27.  CONTINUED 
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THIS  PAGE  IS  BEST  QUALITY  PEA  Cl  lOW 
m>M  CiO>PX  XL  LAG 


RE  ITEI5, 80001  I PLLSUMU  I ,K=3,fc) 
8000  FORMAT!  )X  , A I IPc  l U.  ) ,4X  I I 

1 5 CO  CONTINUE 

GC  TC  25u 

loco  continue 

or  tc  10 

400  continue 

CLOSE! UN IT=5) 

STOP 

END 


SUBRPIT INE  DATE  I L ( J , I ACC , B C , VC , S L C P t C , B I N TC , V CC , S I oDO 1 
DIMENSION  I AC! 12) , 8! 121, V! 12), SLOPE (12  I ,BINT(  12) ,VU( 12) ,SIGUl 121 


DATA 

( IAC (Lit 
L=l, 121/ 

B(L)  , 

VI  LI  , 

S LOPE  I L I 

, e I NT  ( L I , 

> VO  IL  1, 

S1GCU 

l 

•B707  ■ , 

130.9, 

232.3, 

5.491 , 

85.525, 

5.22, 

1.62, 

£ 

•B727  •, 

108. 0, 

205.8, 

5.471, 

83.  186, 

6.55  , 

1.9  , 

3 

•B/3  7 • , 

93.0, 

197.0, 

3.424, 

59.729, 

6.54, 

1.90, 

<* 

1 B 74  7 *, 

195.7, 

238.1  , 

5.  5o  3 , 

106.712, 

6.32  , 

1.94, 

5 

•DC-8  •, 

142.3, 

221  .2, 

4.122, 

118.034, 

4.86, 

1.60, 

6 

•DC-9  •, 

93.  3, 

189.6, 

3.  746  , 

77.763  , 

6.24, 

1.90, 

7 

•CC-  10*  , 

lc.5.3. 

232  .3  , 

4.  712, 

93.864, 

7.00, 

1.9C, 

8 

•L10  111, 

155.  3, 

241.1, 

4.983, 

89.051  , 

6.98  , 

1.90, 

9 

•B707H*, 

145.8, 

232.3, 

4.536, 

93. 593, 

5.  52, 

1.60, 

1 

• DC  - 8H  ’ , 

148.4, 

210.2, 

4.915, 

91.254, 

5.72, 

1.80, 

l 

•PA28  ', 

30.0, 

1 10.0, 

2. 5C8, 

46.475, 

4.00  • 

1.70, 

i 

■ LEAR  • , 

35.6, 

154.0, 

1.586, 

66.441 , 

7.51, 

1.90/ 

IACO-1ACI  J) 

BC  = B I J I 
VC=  V 1 J I 

SLCPEC=SLCPE(J) 
BINTC=BINT( Jl 
VUC  = VC  I J I 
SIGOC=SIGDI J) 
RETURN 
ENC 


SUBROUTINE  NDTRIX.P.C) 

C FOR  X.GE.O,  CALCULATES  AREA  UNDER  NORMAL  0 I S TR I BUT IUN  CURVE 

C FOR  X=X  TO  INFINITY 

C IBM  STATICS  PACKAGE  FSS  PAGE  78 

AX  = ABS I X I 

T=  1.0/1  1.0+0.2316419*AX I 
D=0 . 3989423»EXP<-X*X/2.0) 

P = D*T* (((( l.330274*T-l. 82 1256) *T  + l .781478  I *T-0.3565638) *T 
* +0.3153815) 

RE  TURN 
END 


SUBROUTINE  ERF  UN  IR  . R , S I G ,PE  I 
EX  = ( W + R ) /S  I G 
EY=(R-R)/SIG 
CALL  NDTRI EX.PEX.DEX) 

CALL  NDTRI EY.PEY.DEYI 
IFIEY.LT.O.)  GO  TO  20 
PE  = FEY-FEX 


FIGURE  27.  CONTINUED 
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BSS,QUALinmcnc«« 


IBIS  PACE  IS 

FRW*  om  FUWilbHHJ  10  ^ 


GO  TO  21 

20  PE=1.0-(PEY+PEX) 
21  RETURN 

END 


SUBROUTINE  PwSET 
C0HMCN/PC*1ND/PC>.(51,2I 
DATA  A, B, EMU/12. 5, 5. 5, 18.6/ 
SIG=SCRTI2./3.1<>15926535)*EMU 
TSIG2=SIG*S1G*2. 

B2=e*e 
ACB  = A/B 
E HU  1 = 1 ./EMU 
SUM=C. 

UC  2 1=1.51 
wc=i-i 
WC2  = *C*wC 

EXX=EXP[-hC2/TSlG21 
PCM 1 ,1) =EMUl*EXX 
Q =0 . 5 

1F( 1 .oT.61  GO  TC  3 
E PS  = AC3*SCPT  ( B2-V.C2 ) 

CALL  NDTRIEPS/SIG.U.OUM) 

3 PCM  l,2)=C*EXX 

SUH  = SLM+PCR  ( I ,2  ) 

2 CONTINUE 

OC  4 1=1.51 
PCM  1 ,21  = PCwI  1 ,2)  /Sum 
‘t  CCNT1NUE 

KE  TURN 
END 


FIGURE  27.  CONCLUDED 
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SENSITIVITY  ANALYSIS 


In  the  following  analyses,  the  calculations  were  done  by  posi- 
tioning two  aircraft  (B-747  and  DC-9)  at  a 3 - naut ical -mi  1 e separa- 
tion. In  addition,  the  following  aircraft,  the  DC-9,  was  assumed 
to  be  over  the  OM  as  the  probability  of  a vortex  encounter  is 
largest  at  the  OM  location.  VAS-green  conditions  are  assumed,  and 
all  probabilities  are  referenced  to  the  baseline  value  N.  The 
B-747/DC-9  pair  was  used  as  it  gives  the  largest  probabilities  for 
the  airliners. 

The  key  variable  in  the  conservative  safety  analysis  is  the 
factor  f,  the  fraction  of  the  roll-control  authority  needed  to 
counter  the  swirling  velocities  of  a vortex.  In  Section  4.1.2, 
it  was  argued  that  any  value  between  0.5  and  0.756  could  have  been 
selected,  but  one-half  of  0.756  was  used  to  be  conservative.  It  is 
instructive  to  determine  the  value  of  f at  which  N,  the  baseline 
probability,  is  exceeded  for  a jetliner  following  another  jetliner. 
The  probability  of  a vortex  hazard  exceeds  N in  VAS-green  condi- 
tions for  a DC-9  3 nautical  miles  behind  a B-747  only  with  an  f 
less  than  0.32. 

The  flight  speed  enters  into  the  strength  and  hence  the  decay 
rate  of  the  vortices  shed  by  the  lead  or  vortex  - genera t i ng  aircraft, 
and  parameterizes  the  roll-rate  capability  of  the  following  or 
vortex- encountering  aircraft.  Figure  28  shows  how  the  probability 
of  a vortex  hazard  (in  units  of  N)  varies  with  the  flight  speeds 
for  the  DC-9  following  the  B-747  at  3-nautical-mile  spacings  in  the 
vicinity  of  the  OM.  The  two  curves  were  obtained  by  using  the 
nominal  values  of  all  the  parameters  and  varying  just  the  landing 
speed  of  the  B-747  to  get  the  curve  labeled  B-747,  or  the  landing 
speed  of  the  DC-9  to  get  the  curve  labeled  DC-9.  The  B-747  would 
need  to  be  traveling  at  an  inordinately  high  approach  speed  to 
increase  the  probability  to  N.  The  DC-9  would  need  to  be  flying 
about  10  knots  less  than  the  typical  landing  speed  (close  to  stall 


h.igiit  spun  d't/scc) 


FIGURE  28.  NORMALIZED  PROBABILITY  OF  A VORTEX  ENCOUNTER  FOR  A 
DC-9  THREE  NAUTICAL  MILES  BEHIND  A B-747  DURING  VAS-GREEN 
CONDITIONS  AS  A FUNCTION  OF  THE  FLIGHT  SPEEDS  OF  THE  AIRCRAFT. 
The  solid  dots  indicate  the  nominal  values  of  the  flight  speeds 


speed)  before  the  probability  exceeds  N. 

Figure  29  indicates  the  effect  of  varying  the  initial  descent 
speed  of  the  B-747  vortices.  Low  initial  descent  speeds  for  the 
B-747  are  not  uncommon,  but  as  noted  in  Appendix  B,  low  initial 
descent  speeds  (as  compared  with  theoretical  values)  usually  occur 
during  unstable  atmospheric  conditions  when  vortices  have  shorter 
lives.  The  standard  deviation  in  the  descent  speed  (o^)  does  not 
substantially  affect  the  calculations  as  shown  in  Figure  30.  This 
should  not  be  unexpected  as  the  vertical  deviations  of  the  two 
aircraft  about  the  glide-slope  beam  (o  ) are  much  larger  than  the 
deviation  in  descent  distance  (equation  (24)). 

The  size  of  the  hazard  radius  is  related  to  the  magnitude  of 
the  initial  strength  of  circulation  of  the  vortices.  Figure  31 
shows  the  normalized  probability  of  a vortex  hazard  versus  the 
initial  circulation.  Since  the  value  of  the  initial  circulation 
was  obtained  by  fitting  the  maximum  values  of  T'  for  various  air- 
craft wingspans  (Appendix  A),  a large  increase  in  Fq  seems  highly 
unlikely. 

The  wind  model  derived  in  Appendix  C permits  the  calculation 
of  the  probability  of  a hazardous-vortex  encounter  for  no  VAS 
information  and  for  VAS-green  conditions.  Since  the  model  for 
VAS-green  conditions  was  derived  using  a number  of  assumptions,  it 
is  important  to  examine  the  sensitivity  of  the  results  of  the 
safety  analysis  to  the  choice  of  wind  model.  The  probability  of  a 
hazardous -vortex  encounter  for  3-nautical -mi le  separations  ap- 
proaches zero  for  all  jetliner  combinations  if  f is  greater  than 
0.458  regardless  of  the  wind  model  used.  The  difference  between 
the  no-VAS  and  VAS-green  models  is  important  only  for  the  B-727, 

B-  737  , and  DC.-9  following  the  B-747  , and  for  the  PA-28  and  Learjet 
following  most  of  the  airliners. 

The  results  of  this  brief  sensitivity  analysis  show  that  the 
hazardous - vo rt ex  encounter  model  is  not  too  sensitive  to  changes  in 
single  input  parameters.  Although  only  the  DC  - 9 following  the  B-74 
case  has  been  discussed  herein,  it  represents  the  apparent  "worst 
case"  in  terms  of  safety  for  the  airliners  (see  Table  7). 


0. 


INITIAL  VORTEX -DESCENT  SPEED  (l't /sec) 


FIGURE  29.  NORMALIZED  PROBABILITY  OF  A VORTEX  ENCOUNTER  FOR  A 
DC-9  THREE  NAUTICAL  MILES  BEHIND  A B-747  DURING  VAS-GREEN 
CONDITIONS  AS  A FUNCTION  OF  THE  INITIAL  DESCENT  SPEED  OF  THE 
B-747  VORTICES.  The  solid  dot  indicates  the  nominal  value  of 
the  initial  descent  speed. 
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FIGURE  30.  NORMALIZED  PROBABILITY  OF  A VORTEX  ENCOUNTER  FOR  A 
DC-9  THREE  NAUTICAL  MILES  BEHIND  A B-747  DURING  VAS-GREEN 
CONDITIONS  AS  A FUNCTION  OF  THE  STANDARD  DEVIATION  IN  THE 
INITIAL  DESCENT  SPEED  OF  THE  B-747  VORTICES.  The  solid  dot 
indicates  the  nominal  value  of  the  standard  deviation  in  the 
initial  descent  speed. 


FIGURE  31.  NORMALIZED  PROBABILITY  OF  A VORTEX  ENCOUNTER  FOR  A 
DC-9  THREE  NAUTICAL  MILES  BEHIND  A B-747  DURING  VAS-GREEN 
CONDITIONS  AS  A FUNCTION  OF  THE  INITIAL  CIRCULATION  OR  STRENGTH 
Or  THE  B-747  VORTICES.  The  solid  dot  indicates  the  nominal 
value  of  the  initial  circulation. 
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CONFIRMATION  TESTS 
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F.l  LASER-DOPPLER  VELOCIMETER 

There  are  only  two  practical  methods  for  tracking  vortices 
which  are  between  200  and  1000  feet  above  the  ground:  the 
phototheoJolite  tracking  of  smoke-marked  vortices,  and  the  use  of  an 
LDV.  The  former  technique  ^ould  involve  dedicated  aircraft;  only 
those  aircraft  made  available  to  the  program  could  be  investigated, 
and  tracking  would  cease  when  the  smoke  detrained  (not  necessarily 
at  the  death  of  a vortex).  There  are  technique  variants,  such  as: 
using  cameras  with  telephoto  lenses,  using  markers  other  than  smoke 
(chaff  or  snow  would  permit  vortex-tracking  with  a microwave  radar), 
using  a tall  tower  to  disperse  the  smoke,  etc.  An  LDV,  on  the  other 
hand,  can  be  deployed  at  an  airport  on  a noninterfering  basis,  and 
thus,  yield  data  on  various  aircraft  under  operational  conditions. 

An  LDV  measures  the  Doppler  spectrum  of  laser  radiation  back- 
scattered  by  atmospheric  aerosols.  The  instrument  incorporates 
means  to  transmit  the  laser  radiation  to  the  region  of  interest, 
to  collect  the  radiation  scattered  from  the  aerosols,  and  to 
photomix  on  a photodetector  the  scattered  radiation  and  a portion 
of  the  transmitted  beam.  A di fference- frequency  component,  at  the 
Doppler- sh i ft  frequency,  is  generated  at  the  photo  detector,  and  is 
translatable  into  an  along-opt ic- ax  is  wind - veloc ity  component.  For 
a CO2  laser  operating  at  a wavelength  of  10.6  microns,  a Doppler 
shift  of  97.1  kHz  corresponds  to  a 1-knot  1 ine-of-sight  velocity. 
Details  of  the  operation  of  an  LDV  system  may  be  found  in  Ref.  27- 


Figure  32  shows  the  mobile  FAA/TSC  LDV  housed  in  a 24- foot, 
2-1/2-ton,  detachable-body  truck.  The  system  is  deployed  under  the 
approach  to  runway  14R  at  Chicago's  O'Hare  International  Airport. 
Initially,  the  LDV  was  deployed  under  the  approach  to  runway  27R, 
but  the  lack  of  traffic,  in  particular  the  lack  of  B-747s,  prompted 
moving  the  system  to  14R  during  April  1978.  The  data  collection 
and  relevant  analyses  are  described  in  Volume  II  of  this  report. 

In  the  vortex  mode,  the  LDV  is  used  to  measure  the  velocity 
field  of  the  vortices  by  scanning  the  laser  beam  through  various 
ranges  and  elevation  angles.  Doppler  spectra  are  processed  by  a 
minicomputer  system  in  the  van;  the  largest  Doppler  shift  (and  hence, 
the  largest  velocity)  is  interpreted  as  the  maximum  vortex  tangential 
velocity  (the  velocity  at  the  core  radius,  vy  of  equation  (3)).  The 
vortices  are  tracked  by  following  the  location  of  the  largest 
Doppler  shifts,  one  from  each  side  of  the  vortex  center.  Vortex 
strength  is  calculated  using  the  velocity  field  (the  method  is 
similar  to  that  used  in  the  MAVSS ; see  equation  (51)  in  Appendix  A). 

In  the  wind  mode,  the  LDV  beam  is  focused  at  a range  of  inter- 
est, and  scanned  360  degrees  in  azimuth  while  holding  the  angle 
from  the  zenith  constant  (typically  30  degrees).  If  the  horizontal 
wind  remains  constant  during  the  time  of  a complete  rotation 
(approximately  2 seconds) , the  line-of-sight  Doppler-velocity 
component  will  be  a sine  wave  whose  amplitude  is  a measure  of  the 
horizontal  wind  magnitude,  whose  phase  is  a measure  of  the  wind, 
and  whose  average  value  is  a measure  of  the  vertical  speed.  By 
refocusing  the  LDV  to  a series  of  heights,  a profile  of  the  three 
wind  components  can  be  mapped. 

F . 2 ABBREVIATED  TEST  PLAN 

Since  the  VAS  equipment  is  in  place  at  O'Hare,  it  permits 
correlation  of  high  altitude  (200  to  1000  feet  above  ground  level) 
vortex  behavior  with  VAS  outputs.  It  is  important  to  collect 
enough  runs  for  the  data  to  be  statistically  significant,  and  to 
experience  various  meteorological  conditions. 
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The  ambient  wind  conditions  as  measured  by  the  VAS  are  being 
recorded  for  correlation  purposes.  It  is  expected  that  the  con- 
ditions which  indicate  that  reduced  separations  may  be  used 
(green-light  conditions)  are  the  same  conditions  which  lead  either 
to  the  quick  demise  of  vortices  aloft  or  to  the  rapid  translation  of 
vortices  out  of  the  landing  corridor.  Anomalous  vortex  behavior,  if 
it  occurs,  is  expected  only  when  the  VAS  indicates  vortex-based 
separations  should  be  used  (red-light  conditions) . 

Aircraft  altitude  above  the  LDV  van  is  being  measured  to 
determine  the  initial  altitude  of  the  vortices  and,  hence,  to  allow 
the  measurement  of  the  vortex-descent  distance.  A theodolite 
emplaced  near  the  runway  threshold,  :and  pointed  up  the  glide  slope, 
is  being  used  to  measure  the  azimuth  and  elevation  of  the  aircraft 
as  they  pass  over  the  LDV. 

To  characterize  the  atmosphere,  an  acoustic  sounder  and  a 
pyranograph  have  been  deployed  at  O'Hare.  The  acoustic  sounder 
monitors  the  heights  of  atmosphere  structure  fluctuations  and 
inversion  layers.  The  pyranograph  monitors  the  cloud  cover  by 
measuring  the  solar  and  sky  insolation;  knowing  the  cloud  cover 
allows  an  assessment  of  the  atmospheric  stability. 

All  the  data  are  being  compiled  in  a data  base  for  analysis. 

For  each  flyby  the  following  information  is  obtained:  aircraft 
type,  VAS  tape  and  run  number,  LDV  tape  and  run  number,  aircraft 
lateral  position  and  height  with  respect  to  the  LDV  system,  time 
code,  VAS-averaged  winds  near  runway  threshold,  VAS  status  (red 
or  green  light),  winds  aloft  measured  by  LDV  (when  available), 
inversion  height(s)  (if  present),  pyranograph  output,  general 
weather  conditions  (precipitation,  gustiness,  etc.),  position 
and  age  of  each  vortex  detected  by  the  LDV,  strength  (when  possible) 
of  the  vortices,  and  (when  known)  type  of  operation  (IFR,  VFR,  or 
visual  approach) . 

Studies  will  be  undertaken  to  identify  the  statistically 
significant  correlations  among  the  vortex  behavio’',  ambient 
meteorological  conditions,  and  VAS  outputs.  In  particular,  the 
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measured  vortex  translation,  descent  and  decay  will  be  compared 
with  the  respective  model  adopted  in  Section  4.  If  the  data  indicate 
that  any  of  the  models  must  be  modified,  the  encounter  probabilities 
will  be  recalculated  and  the  consequences  (if  any)  delineated. 

The  winds  aloft  are  also  measured  by  the  LDV  and  compared  with 
the  VAS-measured  winds.  The  comparison  is  required  to  verify  the 
cross-wind  model  developed  in  Appendix  C.  A data  base  is  being 
compiled  consisting  of:  time  code,  VAS-averaged  winds  near  the 
appropriate  runway  threshold,  VAS  status  (red  or  green  light), 
and  averaged  winds  at  various  altitudes  above  the  LDV.  Cross-wind 
models  will  be  derived  from  the  data  for  the  "no-VAS"  and  the  "VAS- 
green"  situations,  and  compared  with  the  models  in  Appendix  C 
(equations  (52)  and  (57),  respectively).  If  the  validity  of  the 
models  in  Appendix  C cannot  be  substantiated,  the  probability 
analysis  will  be  revised  using  the  LDV-obtained  models. 
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